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RESEARCH  SUMMARY 

The  mathematical  programming  formulatton  first 
developed  by  the  Integrated  Resource  Planning  Model 
(IRPM)  (Kirby  and  others  1980)  is  useful  in  identifying 
the  spatial  arrangement  and  timing  of  land  manage- 
ment activities  and  road  construction  and  reconstruc- 
tion projects  that  efficiently  implement  management 
objectives.  In  addition  to  IRPM,  currently  this  formula- 
tion can  also  be  developed  using  the  Integrated 
Resource  Analysis  System  (IRAS)  (Jones  and  others 
1990),  or  the  project  option  within  FOR  PLAN  (Johnson 
and  others  1986). 

The  IRPM-type  formulation  is  a  mixed-integer 
mathematical  program;  that  is,  iDOth  continuous  and 
0,1  integer  variables  are  present.  The  integer  vari- 
ables arise  in  the  transportation  portion  of  the  formula- 
tion. Road  construction  and  reconstruction  projects 


are  developed  for  specific  segments  of  proposed  and 
existing  roads.  These  decision  variables  are  formu- 
lated as  integers  that  can  only  assume  the  values  of 
1  (road  project  is  selected)  or  0  (road  project  is  not 
selected).  Fractional  values  for  road  projects  are  not 
allowed  because  they  represent  the  constmction  of 
only  portions  of  road  segments.  This  leads  to  a 
discontinuous  road  network,  which  is  an  impractical 
solution. 

Unfortunately,  mixed-integer  mathematical  program- 
ming problems  are,  in  general,  difficult  to  solve,  and 
this  has  hampered  use  of  the  IRPM-type  formulation. 
Experience  has  shown  that  branch-and-bound  algo- 
rithms for  solving  these  mixed-integer  formulations, 
such  as  the  algorithm  available  at  National  Computer 
Center  at  Fort  Collins  (NCC-FC),  are  viable  for  solving 
only  relatively  small  models.  For  models  typical  in 
real-world  applications  containing  200  or  more  integer 
variables,  the  branch-and-tx)und  approach  is  generally 
cost-prohibitive,  if  optimal  solutions  can  be  obtained 
at  all. 

This  paper  describes  a  heuristic  procedure  for 
solving  the  IRPM-type  formulation.  The  procedure  is 
an  iterative  process  in  which  models  are  initially 
formulated  and  solved  as  continuous  linear  program- 
ming problems  (road  projects  are  represented  by 
continuous  variables).  Then  logical  decision  rules 
(based  on  the  structure  of  the  IRPM-type  formulation) 
are  applied  to  the  continuous  linear  programming  (LP) 
solution  to  round  some  of  the  fractional  road  projects 
to  0  or  1 .  The  objective  in  these  rounding  decisions  is 
to  preserve  feasibility  and  obtain  solutions  as  close  as 
possible  to  optimum.  These  rounding  decisions  are 
incorporated  into  the  LP  matrix  and  another  LP 
solution  is  made.  Based  on  this  new  solution,  round- 
ing decisions  made  in  previous  iterations  are  reviewed 
to  determine  if  any  changes  should  be  made  and, 
secondly,  additional  fractional  projects  are  rounded  to 
0  or  1 .  These  decisions  are  incorporated  into  the  LP 
matrix,  another  LP  solution  is  obtained,  and  so  on. 
The  process  continues  until  no  fractional  road  projects 
remain  in  the  LP  solution  and  no  changes  in  the 
previous  decisions  can  be  found  that  would  improve 
the  value  of  the  objective  function.  Feasible  mixed- 
integer  solutions  for  full-scale  models  can  generally  be 
obtained  within  seven  to  12  iterations. 
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INTRODUCTION 

Linear  programming  applications  for  timber  harvest  scheduling  and 
resource  allocation  have  been  common  in  forestry  since  the  early  1970's. 
Examples  include  FORPLAN  (Johnson  and  others  1986)  and  earlier  models 
such  as  MUSYC  (Johnson  and  Jones  1979),  Timber  RAM  (Navon  1971),  and 
the  Resource  Allocation  Analysis  system  (USDA  Forest  Service  1975).  These 
applications  have  been  oriented  towards  strategic  planning  where  time 
frames  are  quite  long  and  the  areas  analyzed  quite  large.  They  have  been 
used  to  develop  the  standards  and  guidelines  by  which  forest  land  is  to  be 
managed,  as  well  as  output  targets.  They  do  not,  however,  identify  where 
or  when  to  conduct  projects  to  achieve  these  objectives. 

Tactical  planning  is  concerned  with  the  arrangement  and  timing  of 
projects  that  best  meet  the  management  direction  and  objectives  of  a  strate- 
gic plan.  There  are  often  numerous  considerations  in  making  tactical  timber 
harvest  and  transportation  decisions  for  implementing  long-range  strategic 
plans.  These  decisions  address  the  location  of  the  harvest  units,  silvicul- 
tural  methods,  logging  methods,  and  when  and  how  to  best  access  the  har- 
vest units.  Decisions  about  such  things  affect  the  profitability  of  a  timber 
sale,  the  environmental  effects,  and  the  extent  to  which  multiple-use  objec- 
tives are  met.  But  it  does  not  end  there.  Decisions  about  a  sale  also  have 
implications  regarding  future  management  within  the  general  area  sur- 
rounding the  sale.  The  geographic  placement  of  the  harvest  units  in  the 
current  sale  area  will  affect  where  harvest  units  can  be  located  in  the  future 
because  of  impacts  on  critical  watersheds  and  other  cumulative  environmen- 
tal effects,  wildlife  dispersion  objectives,  and  concerns  about  esthetics. 
Current  activities  may  also  affect  the  financial  aspects  of  selling  the  residual 
timber  in  the  future. 

The  Integrated  Resource  Planning  Model  (IRPM)  (Kirby  1980),  the  Inte- 
grated Resource  Analysis  System  (IRAS)  (Jones  and  others  1990),  and  the 
project  option  within  FORPLAN  (Johnson  and  others  1986)  can  be  used  to 
develop  a  formulation  designed  for  tactical  forest  planning.  Tests  have 
indicated  that  this  formulation  is  quite  useful  in  identifying  land  manage- 
ment and  transportation  projects  that  best  meet  specified  management 
objectives  (Jones  and  others  1986). 

One  problem  that  has  hampered  widespread  use  of  the  IRPM-type  formu- 
lation is  the  lack  of  an  efficient,  easy-to-use  solution  procedure.  Some  of  the 
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decision  variables  in  this  formulation  are  integers  that  can  assume  only 
the  values  of  0  or  1.  Experience  has  shown  that  the  branch- and-bound 
approach  for  solving  mixed-integer  formulations  (such  as  the  Functional 
Mathematical  Programming  System  [Sperry  Corporation  1984],  the  math- 
ematical programming  system  at  the  USDA  National  Computer  Center  at 
Fort  ColUns)  are  viable  for  solving  only  relatively  small  models.  For  models 
tjqjical  in  real-world  applications  containing  200  or  more  integer  variables, 
the  branch-and-bound  approach  has  not  been  practical. 

This  paper  describes  a  heuristic  procedure  for  solving  the  IRPM-type 
mixed-integer  formulation.  It  is  based  on  an  approach  originally  proposed 
by  Weintraub  (1982).  The  procedure  has  been  coded  in  a  program  that 
operates  at  the  USDA  National  Computer  Center  at  Fort  Collins  (NCC-FC). 
This  heuristic  solution  procedure  can  process  models  developed  with  IRPM, 
IRAS,  or  FORPLAN  usiag  the  project  analysis  formulation.  (Readers  inter- 
ested in  more  information  about  this  computer  program  should  contact  the 
authors.) 

The  advantage  afforded  by  this  heuristic  solution  procedure  is  that  fea- 
sible mixed-integer  solutions  can  be  found  for  large,  real-world  problems 
with  much  less  computer  time  than  required  by  a  branch-and-bound  algo- 
rithm. Models  in  the  range  of  2,000  to  4,000  rows,  an  equal  number  of 
columns,  and  200  or  more  integer  variables  can  usually  be  solved  within 
20  minutes  of  central  processing  unit  (CPU)  time  on  the  UNISYS  main- 
frame computer  at  NCC-FC.  Costs  range  from  $50  to  $200,  given  the 
current  rate  structure  at  NCC-FC. 

The  heuristic  procedure  does  not  obtain  optimal  mixed-integer  solutions. 
Tests,  however,  comparing  the  heuristic  procedure  with  a  branch-and- 
bound  algorithm  and  with  other  methods  for  tactical  planning  have  shown 
that  it  provides  good  mixed-integer  solutions,  often  within  10  percent  of  the 
objective  function  value  of  the  true  mixed-integer  optimum. 

The  paper  begins  with  a  presentation  of  the  IRPM-t)T)e  formulation  this 
heuristic  process  is  designed  to  solve.  Next,  the  heuristic  process  is  pre- 
sented, but  largely  in  descriptive  terms.  Much  of  the  detail  regarding  com- 
putations is  presented  in  appendices  A-D.  Last,  results  of  the  heuristic 
process  are  compared  with  results  of  a  branch-and-bound  algorithm  and 
with  other  planning  methods. 


THE  IRPM-TYPE  FORMULATION 

The  IRPM-type  formulation  combines  a  land  allocation  and  scheduling 
model  with  a  cost-minimizing  transportation  network  model.  Land  is 
delineated  into  contiguous  tracts,  hereafter  called  polygons.  These  polygons 
may  represent  potential  harvest  units  (as  illustrated  in  fig.  1),  timber 
stands,  portions  of  stands,  or  other  land  delineations  that  result  in  contigu- 
ous tracts.  Land  management  projects,  called  "resource  projects,"  are  then 
formulated  for  the  polygons.  These  projects  vary  by  the  type  and  timing  of 
management  treatments  to  be  applied.  To  access  these  polygons,  a  road 
network  is  developed  as  illustrated  in  figure  2.  This  network  is  divided 
into  "Hnks,"  which  are  bounded  on  the  ends  by  "nodes."  Road  construction 
projects  are  specified  for  hnks  representing  proposed  new  roads,  and  recon- 
struction projects  are  specified  for  Hnks  representing  existing  roads  requir- 
ing improvement  prior  to  some  specified  future  use.  These  projects  vary  by 
the  standard  to  which  the  road  would  be  constructed  and  the  timing  for 
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Figure  1— Example  of  polygons  delineated  for  an  area. 


3 


Figure  2 — Example  of  a  proposed  road 
network  for  an  area. 
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construction.  Traffic  in  the  model  is  assumed  to  originate  at  the  site  ot  land 
management  activities,  and  is  channeled  through  the  network  to  one  or 
more  "final  demand"  nodes  on  the  extremities  of  the  network.  Round- trip 
impacts  and  costs  are  achieved  by  assuming  traffic  returns  via  the  same 
route.  The  basic  structure  of  the  model  is  presented  below. 


Decision  Variables    1.  Resource  Projects: 

Xiit    =    the  fraction  of  the  resource  project  comprised  of  management 
alternative  I  to  be  implemented  on  polygon  i  in  period  t. 
(0<X/i,<l) 

2.  Road  Projects: 

1     if  the  road  hnk  connecting  adjacent  nodes  a  and  6  is 
to  be  constructed  (or  reconstructed)  to  standard  k  in 
period  t.  For  existing  roads,  a  value  of  1  means 
standard  k  (which  has  no  construction  cost)  is 
available  for  use  starting  in  period  t. 

0  otherwise 

3.  Traffic  Flow  Variables: 

T^ah,t  =     amount  of  traffic  type  y  (for  example,  logging  truck  or 

recreation)  on  road  standard  k  flowing  from  node  a  to  node  6  in 
period  t.  (A  return  trip  in  the  opposite  direction  is  implied.) 

T  bc,t  =     amount  of  traffic  type  >' (for  example,  logging  truck  or 

recreation)  on  road  standard  k  flowing  from  node  b  to  node  c  in 
period  t.  (A  return  trip  in  the  opposite  direction  is  implied.) 

Dbf,t  =    amount  of  traffic  type  y  flowing  from  node  b  and  terminating 
at  final  demand  node  fin  period  t.  (This  is,  in  actuality,  a 
special  type  of  traffic  variable  that  denotes  a  terminating  point 
in  the  road  network.) 


Constraints  1.  Mutually  exclusive  constraints  for  road  construction  projects: 

II    Zabt  <  I  (for  all  a6)  (1) 

k  t 

Only  one  road  project  may  be  selected  for  any  road  link,  ab. 

2.  Traffic  flow  equations: 

Each  traffic  flow  equation  equates  the  volume  of  traffic  type  y  in  period  t 
entering  node  b  to  the  volume  of  traffic  type  y  leaving  that  node  in  period  t. 
In  doing  so  the  traffic  generated  by  the  resource  projects  is  routed  through 
the  road  network  to  the  final  demand  locations.  They  are  formulated  as: 

I  I  Vliu(Xut)+  I   I  (T^'it) 
I    i  aeA  k 

=  11  {Tt^)  +  I  {D^)  (for  all  b,y,t)  (2) 

ceC  k  fzF 
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where 


V 

Vh,iit  =    amount  of  traffic  type  y  loaded  by  resource  project  onto 

node  b  in  period  t  (Vbju  is  zero  for  resource  projects  that 

do  not  load  traffic  onto  node  b) 
A      =    the  set  of  all  nodes  a  from  which  traffic  may  travel  (over  one 

road  hnk)  to  node  6 
C      =    the  set  of  all  nodes  c  to  which  traffic  may  travel  (over  one 

road  link)  from  node  b 
F      =    the  set  of  all  final  demand  nodes  /"connected  to  node  6  by  a 

road  hnk. 

3.  Capacity  equations: 
where 

Ey     =    amount  one  unit  of  traffic  type  y  counts  toward  capacity  on 
^  link  ab 

^ab,g  =    capacity  of  standard  k  (measured  in  a  common  unit  of  traffic) 
on  link  ab  in  time  g. 

Capacity  constraints  serve  two  purposes.  First,  if  link  ab  is  to  carry  traffic 

in  period  t,  the  capacity  equations  ensure  one  road  project  is  selected  for 

link  ab  by  period  t  (that  is,  one  of  the  road  construction  projects  in  one  of 

the  capacity  equations  for  link  ab  must  assume  a  value  of  1).  Second, 

because  each  Z  \i  g  is  restricted  to  the  values  of  0  or  1  and  the  sum  of  these 

variables  defined  for  link  ab  cannot  exceed  1,  capacity  constraints  limit  the 

amount  of  traffic  allowed  in  each  period  g  for  road  standard  k  on  link  ab  to 
/->  k 

^  ab,g- 

4.  Mutually  exclusive  constraints  for  resource  projects: 

X  I  X;if  <  1  (for  all  i)  (4) 

I  t 

Only  one  or  the  equivalent  of  one  resource  project  may  be  selected  for  any 
polygon,!*. 

5.  Typical  optional  relationships  among  resource  projects: 

(i)  Mutually  exclusive  relationships  specify  that  the  equivalent  of  one 
project  may  be  selected  from  a  user-defined  set  S: 

I    Xiu<l  (5) 
il,i,t)£S 

(ii)  Contingent  relationships  require  that  at  least  one  project  from  a 
user-defined  set  S  must  be  selected  if  the  user-specified  project 
Xiit*  is  selected: 

Xiit*<      I  Xiu  (6) 
(l,i,t)zS 

(iii)  Companion  relationships  require  that  if  user-specified  project  X/j^ 
is  selected,  then  some  second  user-specified  project  Xm*  must 
also  be  selected  at  the  same  level,  and  vice  versa: 


(7) 


6 


6.  Resource  Project  to  Road  Construction  Triggers: 

Xiu    <       I.        S     Z*fe^  (8) 
iab,k)eS  g=l 

where  S  is  comprised  of  a  set  of  road  projects,  one  of  which  must  be  selected 
ifXiii  is  selected.  These  triggers  are  used  to  ensure  one  or  more  specific  road 
segments  in  the  vicinity  ofXm  are  constructed  ifXin  is  selected.  They  are 
not  used  to  specify  routes  to  a  final  demand  node,  however.  This  is  accom- 
plished via  the  traffic  flow  equations  (equation  2)  and  capacity  equations 
(equation  3)  which,  respectively,  (1)  channel  all  traffic  originating  with  the 
selected  resource  projects  to  one  or  more  final  demand  nodes,  and  (2)  require 
that  a  road  project  be  selected  for  any  Unk  to  carry  traffic. 

7.  Side  Constraints: 

The  term  "side  constraints"  refers  to  a  class  of  user-defined  management 
constraints.  Side  constraints  include  constraints  placed  on  costs  (such  as 
road  construction  costs,  hauling  costs,  or  total  discounted  costs),  on  physical 
outputs  (such  as  timber  yields,  water  production,  or  sediment  produced),  and 
on  net  revenues.  These  constraints  can  contain  any  of  the  decision  variable 
types  or  any  combination  of  them  and  can  be  formulated  as  upper  limit, 
lower  limit,  or  equality  constraints. 

Objective  Various  objectives  may  be  defined.  A  common  one  is  to  maximize  present 

Function  net  value  (PNV): 

maxPNV=  Z  Z  Z  GiuiXiu)  -  Z   Z  Z  nXtizXt) 
lit  k  ab  t 

-  Z   Z   Z    Z  Iab,t(Tab.t^ 


k  y  ab  t 

z  z  z  z 

k  y  ba  t 


-  z  z  z  z  /J;^,  {tI'Ii)  (9) 


where 


Giii    =    net  value  of  management  alternative  /  on  polygon  i  in  period 
t.  For  timber,  Gm  is  computed  as  the  expected  value  of  the 
delivered  logs,  minus  stump-to-truck  costs  and  all  other 
relevant  on-site  costs,  such  as  sale  preparation  and  admin- 
istration, site  preparation  and  brush  disposal,  and  regen- 
eration. (Neither  hauling  cost  or  road  construction  cost 
should  be  deducted  from  the  delivered  log  value  in  G^ 
^         because  these  costs  are  attached  to  other  decision  variables.) 


^ab.t  =    discounted  cost  of  constructing  link  ab  to  standard  k  in 
^  period  t 

^ab,t   =    discounted  cost  of  a  unit  of  traffic  t5^e  y  traveling  from 
node  a  to  6  on  road  standard  k  in  period  t 


^ba,t   =    discounted  cost  of  a  unit  of  traffic  type  y  traveling  from 
node  6  to  a  on  road  standard  k  in  period  t. 

Other  objective  functions  that  may  be  of  interest  are  to  minimize  total 
discounted  cost  (while  meeting  timber  harvest  targets),  to  maximize  timber 
harvest  over  the  planning  horizon  (while  meeting  environmental  objectives), 
or  to  minimize  environmental  impacts  (while  meeting  timber  harvest  targets). 
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THE  HEURISTIC  INTEGER  PROGRAMMING  (HIP) 
PROCEDURE 


The  heuristic  integer  procedure  (HIP)  described  below  attempts  to  find 
the  mix  of  resource  projects,  road  construction  and  reconstruction  projects, 
and  traffic  routing  that  results  in  the  best  objective  function  value,  while 
satisf)dng  the  constraints  that  have  been  imposed.  Both  maximization  and 
minimization  problems  can  be  solved.  All  road  construction  and  reconstruc- 
tion projects  are  set  to  values  of  0  or  1  by  the  procedure.  The  resource 
projects  are  handled  as  continuous  decision  variables  with  an  upper  bound 
of  1. 

The  HIP  procedure  is  an  iterative  process  that  utilizes  a  continuous  LP 
solution  in  each  iteration.  First,  the  mathematical  model  is  reformulated 
as  a  continuous  LP  problem,  and  additionally,  "deviation  variables"  are 
added  to  the  side  constraints  to  avoid  infeasibihties  that  may  result  from 
HIP  rounding  decisions.  Next,  this  reformulated  problem  is  solved  as  a 
continuous  LP  problem.  The  solution  is  processed  by  the  HIP  procedure, 
essentially  a  set  of  decision  rules  based  on  the  IRPM-type  formulation,  to 
determine  which  road  projects  to  round  to  values  of  0  or  1.  These  decisions 
are  incorporated  into  the  LP  matrix,  and  another  continuous  LP  solution  is 
obtained.  The  new  solution  is  processed,  and  the  resulting  decisions  are 
incorporated  into  the  LP  matrix,  and  so  on.  The  process  continues  until  all 
the  road  projects  have  been  assigned  a  value  of  0  or  1,  and  no  further  modi- 
fications can  be  found  that  improve  the  objective  function  value  without 
violating  one  or  more  constraints.  This  generally  requires  from  seven  to 
12  iterations.  The  matrix  modifications  and  decision  rules  are  described 
below. 

Formulation  Several  modifications  are  made  to  the  mathematical  formulation 

Modifications  described  earlier  to  prepare  a  Hnear  programming  matrix  for  the  HIP 

solution  procedure.  First,  the  road  construction  and  reconstruction  projects 
are  defined  as  continuous  decision  variables  that  can  assume  any  nonneg- 
ative  value.  This  is  necessary  to  permit  solving  a  model  as  a  linear  pro- 
gramming problem. 

Second,  a  set  of  "deviation  variables"  are  added  to  each  side  constraint. 
These  deviation  variables  prevent  the  continuous  LP  from  becoming 
"infeasible"  after  variables  have  been  rounded  to  values  of  0  and  1  by  the 
HIP  program.  (When  an  LP  problem  is  determined  to  be  infeasible,  most 
LP  solvers  report  only  the  information  for  the  LP  iteration  in  which  the 
infeasibility  was  detected.  Because  this  is  normally  only  one  of  many 
possible  infeasible  solutions,  these  solutions  do  not  contain  reliable  infor- 
mation for  making  heuristic  decisions.)  The  formulation  for  an  upper 
bound  (less  than  or  equal  to)  constraint  is: 

[Original  Maximization  Objective  Function]  -  C,x  (Dn)  -  Cih  {^ih) 

[Original  Left-Hand  Side  of  Constraint  i]        ADu)        AD  in)  ^  RHSi 

(DiL)  ^  BiL 

where 

DiL    =  lower-cost  deviation  variable  for  constraint  i 
Djjj    =  higher-cost  deviation  variable  for  constraint  i  ■ 
RHSi  =  right-hand-side  value  for  constraint  i 
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CiL  =  penalty  assigned  to  Di£, 
Ciff  =  penalty  assigned  to  Duj 
BiL    =  upper  bound  placed  on  Dj£^ 

To  illustrate,  if  the  original  left-hand  side  of  a  constraint  were  to  exceed 
RHSi,  DiL  and/or  Djjj  would  assume  a  nonzero  value,  thus  avoiding  an  in- 
feasibility  in  the  LP  solution.  The  same  formulation  is  used  for  minimi- 
zation objective  functions,  with  the  exception  that  the  objective  function 
penalties,  Cn  and  CiH^  carry  positive  rather  than  negative  signs. 

The  penalty  for  violating  a  constraint  is  imposed  in  a  step- wise  manner, 
that  is,  CiH  >  CiL-  The  logic  is  that  some  amount  of  constraint  violation  can 
usually  be  tolerated  because  of  the  inaccuracies  and  random  nature  of  some 
constraint  coefficients.  The  amoimt  that  can  be  tolerated,  Bn,  is  specified  by 
the  user  (the  default  is  10  percent  oiRHSi).  Any  constraint  violation  above 
this  amoimt  is  forced  into  the  higher-cost  deviation  variable,  Djj,  whose 
coefficient,  Cj^,  is  1.2  times  Cil.  When  this  occurs,  the  heuristic  considers 
the  problem  to  be  infeasible.  The  heiuistic  solution  output  reports  the 
objective  function  value  both  with  and  without  the  deviation  variables. 

The  lower-cost  penalty,  C,£„  is  calculated  as  follows: 

Cii  =1.5  times  the  larger  of:  Rmax  or  1 
where 

Rmax  =  the  maximum  I  Cj/aij  I  (for  all  j  where  Oy  ^  0) . 

Here,  Cj  is  the  objective  function  coefficient  for  the  Jth.  decision  variable,  and 
Oy  is  the  constraint  coefficient  for  the  j'th  variable  in  constraint  i.  ex- 
ceeds the  contribution  per  imit  of  constraint  i  that  each  decision  variable 
makes  in  the  objective  function,  so  that  deviation  variables  are  in  solution  at 
a  nonzero  value  only  when  necessary  to  maintain  feasibility  in  a  continuous 
LP  solution.  In  the  absence  of  interdependence  among  decision  variables 
this  would  be  accomplished  by  any  Cn  value  exceeding  Rmax-  Relationships 
such  as  companion  relationships,  however,  may  cause  c^/Oy  to  underesti- 
mate the  objective  function  impact  per  unit  of  side  constraint  i.  Rmax  is 
multiphed  by  1.5  to  allow  for  such  relationships. 

If  there  are  no  decision  variables  for  which  both  cj  and  Cy  are  nonzero, 
then  CiL  is  calculated  as: 

C,x  =1.5  times  the  larger  of:  Smax  oi*  1 

where 

Smax  -  the  maximum  I  cj  I  (over  all  f). 

Low-  and  high-cost  deviation  variables  are  also  added  to  lower  bound 
(greater  than  or  equal  to)  constraints.  The  only  difference  is  that  positive 
signs  are  placed  on  the  constraint  coefficients  for  the  deviation  variables. 
Four  deviation  variables  are  added  for  equality  constraints,  a  pair  of  lower 
and  higher  cost  variables  for  exceeding  the  RHS  value  and  a  pair  for  under- 
achieving the  RHS  value.  Four  deviation  variables  are  also  added  when  a 
side  constraint  contains  both  an  upper  and  lower  bound  (called  ranges). 

After  the  formulation  modifications  have  been  made,  the  first  continuous 
hnear  programming  solution  is  obtained. 


Iteration  1:  Actions      Iteration  1  is  a  capacity-modification  iteration.  The  coefficients  associated 
Following  the  First   with  the  road  projects  in  the  capacity  equations  (equation  3  in  the  general 
LP  Solution  formulation)  are  lowered  to  create  larger  firactional  values  for  the  road  projects 
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this  ratio  is  small.  As  a  result,  something  less  than  the  full  cost  of  con- 
structing link  ab  is  accounted  for  in  the  continuous  LP  solution.  This  is 
often  the  case  in  the  first  LP  solution  because  initially  (as  discussed  earher) 
these  coefficients  represent  the  maximiun  traffic  capacity. 

The  capacity  coefficients  are  reduced  to  more  closely  approximate  the 
observed  traffic  volume  on  the  network.  As  a  result,  road  construction 
projects  enter  into  subsequent  LP  solutions  at  larger  fractional  values,  so 
that  more  of  the  road  construction  costs  are  included.  This  provides  im- 
proved LP  solutions  on  which  to  base  future  heuristic  decisions.  Once  a 
decision  has  been  made  to  round  a  road  construction  project  for  link  a6  to  1, 
the  original  capacity  coefficients  for  that  link  are  restored  to  Hmit  traffic  to 
the  intended  capacity.  (See  appendix  A  for  a  detailed  presentation  of  the 
capacity  adjustment  computations.) 

Other  LP  Matrix  Modifications — ^Two  additional  changes  are  made  to 
the  LP  matrix  prior  to  obtaining  the  next  LP  solution.  First,  the  constraints 
that  set  the  sum  of  the  road  projects  for  a  hnk  to  less  than  or  equal  to  1 
(equation  1)  are  made  into  free  (nonconstraining)  rows.  This  permits  road 
construction  projects  to  assume  values  greater  than  1,  a  change  made  nec- 
essary by  the  capacity  adjustment  discussed  above.  When  the  capacity 
coefficients  are  reduced  to  increase  the  fractional  value  of  road  projects  in 
solution,  there  is  a  chance  that  the  volume  of  traffic  to  be  routed  over  a  link 
in  a  later  solution  will  exceed  the  new  (smaller)  capacity.  Freeing  these 
constraints  allows  this  to  occur. 

Second,  the  objective  fimction  coefficients  for  the  traffic  variables  are 
divided  by  100.  The  rationale  for  this  modification  is:  In  the  continuous  LP 
solutions,  two  types  of  transportation  costs  are  considered  in  routing  traffic: 

(1)  the  cost  associated  with  the  flow  variables  (haul  cost),  and  (2)  the  cost 
associated  with  road  construction  or  reconstruction.  The  continuous  LP 
formulation  does  a  very  good  job  of  routing  traffic  so  as  to  minimize  haul 
cost,  because  the  traffic  variables  were  originally  formulated  as  continuous 
decision  variables.  But  for  reasons  discussed  under  capacity  adjustment, 
road  construction  costs  may  not  be  included  in  their  entirety  (because  the 
road  project  is  in  the  solution  at  a  fractional  value).  As  a  result,  in  the 
continuous  LP  solutions  more  weight  is  placed  on  haul  costs  than  road 
construction  costs  in  routing  traffic.  This  is  offset  in  the  following  LP 
solution  by  reducing  the  objective  function  coefficients  for  the  flow  variables. 
After  the  second  LP  solution  is  obtained,  the  original  objective  function 
coefficients  are  restored.  But  the  information  about  the  least-road  con- 
struction cost  routes  is  captured  via  another  capacity  modification  in  the 
next  iteration. 

The  primary  action  taken  in  this  iteration  is  a  second  capacity  adjustment. 
The  purpose  of  this  adjustment  is  to  further  modify  capacity  on  the  Hnks  to 
better  approximate  the  actual  traffic  volumes,  while  at  the  same  time  en- 
couraging future  traffic  over  hnks  carrying  relatively  large  traffic  volumes  in 
the  current  LP  solution. 

Adjustment  of  Capacity  Coefficients — ^The  capacity  for  Hnk  ab  is  modified 
to  the  larger  of:  (1)  the  average  observed  traffic  flow  on  the  network,  or 

(2)  the  observed  traffic  flow  on  link  ab.  This  adjustment  encourages  traffic  ia 
subsequent  LP  solutions  to  travel  over  the  Hnks  carrying  relatively  high  vol- 
imies  of  traffic  in  the  current  LP  solution  where  the  effects  of  road  construction 


Iteration  2:  Actions 
Following  the 
Second  LP  Solution 
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in  the  subsequent  continuous  LP  solutions.  These  larger  values  for  road 
projects  provide  an  improved  basis  from  which  to  make  rounding  decisions. 
After  a  road  project  has  been  rounded  to  either  0  or  1,  the  original  coeflBdents 
are  reinstated.  In  addition,  the  first  LP  solution  is  analyzed  in  this  iteration 
to  see  if  processing  should  stop  because  the  mixed-integer  problem  is  either 
infeasible,  unbounded,  or  that  the  optimal  solution  is  to  do  nothing. 

Checks  for  Infeasible,  Unbounded,  or  Null  Solutions — If  the  first 
continuous  LP  solution  is  infeasible,  processing  is  stopped  because  the 
mixed-integer  problem,  which  is  more  restrictive,  must  also  be  infeasible. 
There  are  two  checks  for  infeasibihty.  First,  the  LP  solution  status  is 
checked  to  determine  if  it  is  "infeasible."  Second,  the  LP  solution  is  checked 
to  determine  if  a  higher-cost  deviation  variable  is  in  solution  with  a  nonzero 
value.  This  occm^  when  the  LP  solver  does  not  find  a  continuous  solution 
within  the  bounds  placed  on  the  side  constraints. 

If  the  first  continuous  LP  solution  is  unbounded,  processing  is  stopped 
because  this  indicates  an  error  in  the  LP  matrix.  All  the  decision  variables 
have  upper  limits  imposed  by  the  relationships  in  the  IRPM-type  formu- 
lation. The  road  construction  projects  are  hmited  to  the  equivalent  of  one 
per  road  link  (equation  1);  the  resource  projects  are  limited  to  the  equi- 
valent of  one  per  polygon  (equation  4);  and  the  traffic  variables  cannot 
exceed  the  capacity  coefficients  in  equation  3.  Thus,  unbounded  LP 
solutions  can  be  obtained  only  if  there  is  something  wrong  with  the  LP 
matrix. 

If  no  projects  were  selected  in  the  first  continuous  LP  solution,  processing 
is  stopped.  In  this  case  the  optimal  mixed-integer  solution  is  to  do  nothing. 

Adjustment  of  Capacity  Coefficients — If  the  linear  programming 
solution  is  feasible,  optimal,  and  nonnuU,  then  the  capacity  coefficients  in 
the  LP  matrix  are  adjusted  downward  to  match  the  largest  observed  flow 
for  each  road  standard. 

The  capacity  coefiBdents,  C      are  assodated  with  the  road  construction 
projects  in  the  capadty  equations  (equation  3).  To  understand  the  logic  behind 
adjusting  them,  one  must  first  imderstand  the  role  they  play.  Initially,  these 
coefiBdents  represent  the  maximum  trafific  volume  permitted  on  Hnk  ab  if 
constructed  to  road  standard  k.  Given  the  volume  of  trafific  normally  found 
on  forest  roads,  this  limit  is  t3T)ically  much  larger  than  the  actual  volume  of 
trafific. 

These  coefiBdents  play  another  role  in  the  continuous  LP.  They  determine 
the  amount  of  road  construction,  and  therefore  the  road  construction  cost  that 
is  induded  in  a  continuous  LP  solution.  This  is  most  easily  explained  via  a 
simplified  example.  The  capadty  equation  for  Hnk  ab,  having  only  one  timing 
option  and  one  road  standard  option,  would  be  written  as: 

1.0  (Tabl)  -  CAP  (abl)  <  0 

where  Tabl  is  the  traffic  volume  over  link  ab,  abl  is  the  road  construction 
project  (formulated  as  a  continuous  variable),  and  CAP  is  the  capacity  co- 
efifident.  Solving  for  abl  yields: 

abl  >  {\.QICAP)Tabl 

Since  abl  imposes  a  cost  in  both  cost  minimization  and  profit  maximization 
objective  functions,  the  optimization  process  in  the  continuous  LP  will 
assign  a  value  to  abl  as  small  as  possible.  This  value  will  be  Tabl  I  CAP. 
When  capacity  is  substantially  larger  than  the  actual  quantity  of  trafific. 
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costs  are  magnified  relative  to  haul  costs.  (See  appendix  B  for  a  detailed 
presentation  of  the  capacity  adjustment  computations.) 

Other  LP  Matrix  Modifications — One  additional  change  is  made  to  the 
LP  matrix  prior  to  obtaining  the  next  LP  solution.  The  original  objective 
function  coefficients  for  the  traffic  decision  variables  are  reinstated  in  the 
LP  matrix. 

Iteration  3:  Actions  Iteration  3  is  the  first  of  the  iterations  in  which  rounding  decisions  are  made. 
Following  the  Rounding  Procedure:  Round  Selected  Road  Projects  to  1— The  first 

Third  LP  Solution      step  in  the  rounding  procedure  is  to  develop  a  fist  of  candidate  road  projects 

for  rounding.  Each  Hnk  having  flow  in  the  current  continuous  LP  solution 
provides  one  candidate  project.  The  project  selected  has  (1)  the  lowest  road 
standard  that  will  handle  the  volume  of  flow  over  all  periods  on  the  Hnk 
and  (2)  a  period  of  implementation  that  corresponds  to  the  first  period  in 
which  there  is  a  significant  amount  of  traffic  on  that  link  (flow  on  the  link 
exceeds  5  percent  of  the  average  flow  over  all  finks  in  that  period).  In 
addition,  the  road  project  selected  must  be  compatible  with  any  resource 
projects  that  may  have  been  fixed  to  1  by  the  user.  For  example,  if  a  fixed 
resource  project  would  require  (through  a  resource  project-to-road  trigger, 
equation  8)  that  link  ab  be  constructed  in  period  1,  but  no  flow  occurs  until 
period  2,  the  period  1  road  project  would  become  the  candidate. 

Next,  an  index  for  ranking  candidate  road  projects  is  calculated  for  each 
candidate  project.  Projects  are  ranked  in  ascending  order  by  index  value. 
(Projects  with  smaller  index  values  are  preferable  to  projects  with  higher 
index  values.)  The  index  is  comprised  ofthesimiofthe  following  components: 

1.  The  contribution  of  the  road  project  to  the  side  constraints  relative  to 
the  slack  (unused  space)  in  the  side  constraints  weighted  by  1  minus  the 
activity  level  (solution  value)  for  the  road  project  in  the  current  LP  solution. 
The  greater  this  component,  the  larger  (less  desirable)  the  index.  The 
weights  cause  the  contribution  to  the  index  to  be  less,  the  larger  the  activity 
level.  Thus,  the  important  road  projects  in  the  continuous  LP  (as  measured  by 
the  magnitude  of  the  activity  level)  have  less  penalty  for  side  constraint  impact 
than  do  road  projects  in  solution  at  lower  fractional  values. 

2.  The  amount  of  flow  on  link  ab  relative  to  the  average  flow  on  the 
network  in  the  current  LP  solution,  multipHed  by  -1.  The  more  negative 
this  component,  the  more  desirable  the  index.  This  component  uses  quantity 
of  traffic  to  impute  a  measure  of  importance  for  the  links. 

3.  The  relative  cost  of  the  road  project,  multiplied  by  a  weight  of  1  minus 
the  activity  level  for  the  road  project  in  the  current  LP  solution.  The 
greater  this  component  the  less  desirable  the  index.  As  in  1  above,  the 
weight  in  this  component  reduces  the  amount  of  penalty  imposed  in  the 
index  for  projects  having  larger  activity  levels.  Thus,  the  more  important 
projects  (as  measured  by  the  size  of  their  activity  level)  have  less  penalty 
for  construction  cost. 

Next,  determine  which  candidate  projects  connect  to  the  "developed 
portion  of  the  road  network."  The  "developed  portion  of  the  road  network" 
consists  of  road  projects  representing  existing  roads,  plus  projects  repre- 
senting road  construction  or  reconstruction  that  were  previously  rounded 
to  1.  Only  road  projects  connecting  to  the  "developed  portion  of  the  road 
network"  may  be  rounded  to  1.  This  gives  rise  to  a  continuous  developed 
network  that  connects  to  the  demand  nodes. 
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Last,  from  the  set  of  road  projects  connecting  to  the  "developed  portion  of 
the  road  network,"  round  to  1  the  project  with  the  best  (lowest)  index,  pro- 
vided that  project  is  feasible  with  each  of  the  side  constraints.  If  the  project 
is  not  feasible,  select  the  connecting  project  with  the  next  best  (lowest)  in- 
dex, providing  it  is  feasible.  At  most,  no  more  than  50  percent  plus  one  of 
the  original  candidate  projects  can  be  rounded  to  1  in  iteration  3  or  4.  This 
is  increased  to  75  percent  plus  1  beginning  with  iteration  5.  These  cut-off 
points  are  a  compromise  between  making  as  many  decisions  as  possible  in 
an  iteration  to  reduce  the  total  number  of  iterations  and  avoiding  making 
too  many  decisions  which,  as  experience  has  shown,  can  result  in  poor 
mixed-integer  solutions.  (See  appendix  C  for  a  detailed  presentation  of 
the  rounding  procediu^e.) 

The  feasibility  of  rounding  a  candidate  road  project  with  regard  to  side 
constraint  i  is  checked  by  adding  the  impact  of  that  road  project  to  the 
cumulative  impact  from  projects  previously  rounded  in  the  current  iteration, 
and  comparing  that  total  to  the  "allotted  space"  in  side  constraint  i.  The 
"allotted  space"  is  total  amount  of  impact  from  rounding  projects  that  is 
permitted  in  the  current  iteration.  If  the  "allotted  space"  is  not  exceeded, 
the  project  is  feasible. 

Calculation  of  "allotted  space"  begins  with  the  "actual  space,"  which  is 
an  estimate  of  the  actual  space  available  in  a  side  constraint  for  absorbing 
impacts  from  rounding  fractional  road  projects  to  1.  The  "actual  space"  is 
calculated  immediately  following  an  LP  solution.  First,  the  total  impact 
from  the  previously  rounded  road  projects,  the  traffic  variables,  and  the 
resource  projects  is  calculated  by  subtracting  the  impacts  of  the  fractional 
road  projects  from  the  activity  level  for  side  constraint  i  in  the  LP  solution. 
Next,  this  total  impact  is  subtracted  from  the  right-hand-side  value  of  the 
constraint  to  compute  the  estimated  "actual  space." 

The  "allotted  space"  is  then  computed  by  adjusting  "actual  space"  by  an 
adjustment  percentage.  This  percentage  is  composed  of  two  factors  multi- 
plied together.  The  first  factor  varies  the  direction  and  amount  of  adjust- 
ment based  on  the  relative  importance  of  road  construction  in  side  con- 
straint i.  When  road  construction  comprises  the  majority  of  the  impacts  in 
side  constraint  i,  this  factor  would  reduce  "actual  space"  by  approximately 
50  percent.  The  amount  of  reduction  becomes  less  as  the  relative  impor- 
tance of  road  construction  in  side  constraint  decreases,  reaching  zero  when 
roads  comprise  10  percent  of  the  side  constraint  impact.  These  reductions 
are  designed  to  avoid  committing  too  much  of  the  available  space  in  a  side 
constraint  in  any  one  iteration  when  roads  play  an  important  role  in  the 
side  constraint.  (Experience  has  shown  that  objective  function  value  of  the 
resulting  HIP  solution  suffers  when  a  large  portion  of  the  rounding  decisions 
are  made  in  any  one  iteration.)  When  roads  comprise  less  than  10  percent 
of  the  side  constraint  impact,  this  factor  applies  a  percentage  increase  to 
"actual  space."  The  maximimi  upward  adjustment  of  10  percent  is  reached 
when  roads  comprise  less  than  approximately  8  percent  of  the  side-constraint 
impact.  In  these  cases,  roads  play  such  a  minor  role  that  impacts  exceeding 
"actual  space"  can  be  accommodated  in  subsequent  LP  solutions  by  adjust- 
ments in  the  resource  projects  and  traffic  flow  variables  selected. 

The  second  factor  phases  out  the  first  factor  as  the  number  of  iterations 
increase  so  that  the  "actual  space"  in  the  side  constraints  eventually  be- 
comes available.  The  reduction  of  the  first  factor  begins  in  iteration  4,  the 
second  rounding  iteration,  and  increases  in  subsequent  iterations  until  it  is 
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completely  phased  out  in  iteration  8.  Thus,  beginning  with  iteration  8, 
"allotted  space"  equals  "actual  space."  (See  appendix  D  for  a  detailed 
presentation  of  the  feasibiUty  computations.) 

Other  LP  Matrix  Modifications — Two  additional  changes  are  made 
to  the  LP  matrix  prior  to  obtaining  the  next  LP  solution.  First,  the  original 
capacity  coefficients  are  reinstated  in  the  matrix  for  those  links  having  a 
road  construction  project  rounded  to  1  this  iteration.  Second,  the  mutually 
exclusive  constraints  (equation  1)  associated  with  the  links  having  a  project 
rounded  to  1  are  reinstated  as  less  than  constraints  (they  were  changed  to 
nonconstraining  rows  in  the  first  iteration). 

Iteration  4:  Actions       In  iteration  4  and  all  subsequent  iterations,  the  traffic  volumes  on  the 
Following  the  links  that  contain  firactional  road  projects  are  compared  with  the  current 

Fourth  LP  Solution    capacity  coefficients  on  those  links.  If  the  difference  exceeds  a  threshold 

value,  this  iteration  becomes  exclusively  a  capacity  adjustment  iteration. 
If  the  threshold  is  not  exceeded,  this  iteration  reviews  the  rounding  deci- 
sions made  in  previous  iterations,  and  then  proceeds  to  the  procedure  for 
rounding  fractionsil  road  projects. 

Adjust  Capacity  Coefficients  if  Needed — Beginning  with  iteration  4, 
the  observed  flow  on  the  links  having  fractional  road  projects  is  compared 
with  the  current  capacity  on  those  hnks.  If  the  observed  flows  are,  on  the 
average,  less  than  10  percent  of  capacity  for  the  standard  that  matches  the 
observed  quantity  of  traffic  (the  lowest  road  standard  that  can  accommo- 
date the  flow  on  those  links  over  all  periods)  and  capacity  was  not  adjusted 
in  the  previous  iteration,  a  capacity  adjustment  is  made.  This  adjustment 
is  identical  to  the  iteration  2  procedure  and  is  described  in  detail  in  appendix  B. 

It  may  not  be  obvious  why  further  capacity  adjustments  are  necessary. 
The  capacity  adjustment  in  iteration  2  is  based  on  the  average  traffic  flow 
over  links  having  fractional  road  projects  at  that  time.  As  the  rounding 
iterations  proceed,  the  road  projects  rounded  first  are  those  closest  to  the 
demand  nodes  carrying  the  larger  quantities  of  traffic.  The  remaining  links 
having  fractional  road  projects  tend  to  be  in  the  smaller  tributaries  of  the 
network  and  carry  a  smaller  amount  of  traffic.  If  the  flow  on  these  Hnks  is 
small  relative  to  capacity,  it  is  advantageous  to  adjust  capacity  downward 
for  the  reasons  described  earlier — ^larger  fractional  values  for  road  projects 
in  later  iterations,  which  provide  a  better  basis  for  roimding  projects. 

If  capacity  for  the  fi*actional  links  is  adjusted,  the  modified  capacity 
coefficients  are  incorporated  into  the  LP  matrix,  and  another  LP  solution 
is  obtained  for  the  next  iteration.  Alternatively,  if  the  observed  flows  do 
exceed  10  percent  of  capacity  on  the  average,  no  capacity  adjustment  is 
undertaken  and  the  iteration  proceeds  as  described  below. 

Change  Road  Standard  Selections? — Road  projects  rounded  to  1  in 
previous  iterations  are  reviewed  to  determine  if  they  should  be  replaced 
with  projects  having  either  a  lower  or  higher  road  standard.  The  criteria 
are  as  follows: 

1.  Lower  the  road  standard  for  Unk  ab  if  the  largest  traffic  flow  in  any 
one  period  is  less  than  the  capacity  of  a  lower  road  standard  in  that  period. 
Tentatively  select  the  lowest  standard  whose  capacity  exceeds  the  observed 
traffic  flow  on  link  ab  in  each  period.  If  the  lower-standard  project  is  fea- 
sible with  regard  to  the  side  constraints,  round  it  to  1  and  fix  to  0  all  other 
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projects  defined  for  that  link.  If  the  tentatively  selected  project  is  infea- 
sible,  do  not  make  any  change  in  standard. 

2.  Raise  the  road  standard  selection  for  link  ab  if  the  traffic  flow  in  any 
period  equals  the  capacity  for  the  road  project  previously  rounded  to  1. 
Tentatively  select  the  next  higher  standard.  If  that  project  is  feasible  with 
regard  to  the  side  constraints,  round  it  to  1  and  fix  to  0  all  other  projects 
defined  for  that  Unk.  If  the  tentatively  selected  project  is  infeasible,  do  not 
make  any  change  in  standard. 

Rounding  Procedure:  Round  Selected  Road  Projects  to  1 — If  modifi- 
cations were  made  to  "the  developed  portion  of  the  road  network"  during  this 
iteration  (for  example,  change  in  road  standard),  and  those  modifications 
resulted  in  a  critical  change,  then  this  iteration  stops  without  rounding  any 
additional  fractional  projects.  A  change  is  considered  critical  if: 

1.  Any  side  constraint  that  was  within  10  percent  of  its  bounds  before 
modifying  "the  developed  portion  of  the  road  network"  is  no  longer  within  10 
percent  of  its  boimds  after  the  modifications. 

2.  Any  side  constraint  that  was  not  within  10  percent  of  its  bounds  before 
modifying  "the  developed  portion  of  the  road  network"  is  now  within  10 
percent  of  its  bounds  after  the  modifications. 

(See  appendix  D  for  a  detailed  presentation  of  this  feasibility  testing.)  Modi- 
fications in  "the  developed  portion  of  the  network"  resulting  in  critical 
changes  can  have  a  significant  effect  on  the  selection  of  resoiu-ce  projects, 
firactional  road  projects,  and  traffic  routing.  So,  if  either  of  these  conditions 
are  present,  the  LP  is  solved  again  and  rounding  decisions  are  postponed. 

If  modifications  in  the  road  standard  selections  do  not  produce  a  critical 
change,  the  rounding  procedure  described  earher  for  iteration  3  is 
undertaken. 

Other  LP  Matrix  Modifications — These  are  the  same  as  the  changes 
made  in  iteration  3. 

Actions  taken  in  iteration  5  are  identical  to  iteration  4  with  the  exception 
of  an  additional  check  to  determine  if  the  timing  for  construction  should  be 
postponed  to  a  later  period. 

Adjust  Capacity  Coefficients  if  Needed — ^This  step  is  the  same  as  in 
iteration  4. 

Change  Road  Standard  Selections? — This  step  is  the  same  as  in 
iteration  4. 

Postpone  Road  Construction? — Road  projects  rounded  to  1  in  previous 
iterations  are  reviewed  to  determine  if  they  should  be  replaced  with  projects 
implementing  road  construction  in  a  later  period.  Postponing  construction 
is  attempted  if  in  the  period  of  construction  the  amount  of  traffic  on  link  ab 
is  less  than  5  percent  of  the  average  flow  on  the  links  in  the  network.  Tenta- 
tively select  the  project  of  the  same  standard  whose  construction  period 
coincides  with  the  earUest  period  in  which  flow  exceeds  5  percent  of  the 
average  network  flow.  If  this  project  is  feasible  with  regard  to  the  side 
constraints,  round  it  to  1,  and  set  aU  other  projects  for  the  link  to  0.  If  this 
project  is  infeasible  with  regard  to  the  side  constraints,  make  no  change. 
If  there  is  no  later  period  in  which  flow  exceeds  5  percent  of  the  average 


Iteration  5:  Actions 
FoUowing  the  Fifth 
LP  Solution 
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network  flow,  no  change  is  made  in  the  postponing  routine.  This  situation, 
however,  is  addressed  in  the  "close  links"  routine  beginning  in  iteration  6. 

Rounding  Procedure:  Round  Selected  Road  Projects  to  1 — This 
step  is  the  same  as  in  iteration  4. 

Other  LP  Matrix  Modifications — ^These  are  the  same  as  the  changes 
made  in  iteration  3. 


Iterations  6+: 
Actions  Following 
the  Sixth  to  Last 
LP  Solution 


This  iteration  is  identical  to  iteration  5  with  the  exception  of  an  additional 
check  to  determine  if  any  link  previously  selected  for  road  construction 
should  be  closed. 

Adjust  Capacity  Coefficients  if  Needed — This  step  is  the  same  as  in 
iteration  4. 

Change  Road  Standard  Selections? — This  step  is  the  same  as  in 
iteration  4. 

Postpone  Road  Construction? — This  step  is  the  same  as  in  iteration  5. 

Close  Links? — Road  projects  rounded  to  1  in  previous  iterations  are 
reviewed  to  determine  if  sufficient  traffic  is  being  carried  to  warrant  the 
road  project.  If  the  amount  of  traffic  on  link  ab  in  each  period  is  less  than 
5  percent  of  the  average  flow  on  the  links  in  the  network,  that  link  is  closed. 
This  is  accomplished  by  setting  all  road  projects  pertaining  to  that  link  to 
zero  in  the  LP  matrix. 

Rounding  Procedure:  Roimd  Selected  Road  Projects  to  1 — This 
step  is  done  the  same  as  in  iteration  4. 

Other  LP  Matrix  Modifications — These  changes  are  the  same  as  those 
made  in  iteration  3. 

Stopping  Rule — The  solution  process  stops  when  the  current  LP  solu- 
tion contains  integer  values  for  all  road  projects  and  there  are  no  desirable 
changes  to  be  made  to  previous  rounding  decisions  (changing  standard, 
postponing,  or  closing  hnks).  This  generally  occurs  within  seven  to  12 
iterations.  The  final  LP  solution  is  the  integer  solution  for  roads  having 
the  best  objective  function  value,  while  satisfying  the  side  constraints.  In 
nearly  all  cases  this  is  the  last  LP  solution  made  in  the  heuristic  process. 

It  is  possible  that  the  heuristic  solution  process  may  not  be  able  to  find  a 
feasible  mixed-integer  solution  (all  higher-cost  deviation  variables  equal  0) 
even  though  a  continuous  optimal  solution  was  found  in  the  first  iteration. 
This  could  be  because  there  is  no  feasible  mixed-integer  solution  (mixed- 
integer  solutions  are  more  restrictive  than  continuous  solutions).  Or  it 
may  be  that  a  mixed-integer  optimal  exists  but  was  simply  not  found  by 
the  heuristic  procedure.  In  any  case,  when  this  occurs  the  solution  process 
reports  the  best  mixed-integer  solution  it  can  find.  The  value  assigned  to 
the  deviation  variables  measures  the  extent  to  which  side  constraints  are 
violated. 


Comparison  With 
Optimal  Mixed- 
Integer  Solutions 


SOLUTION  VALIDATION 

One  method  to  validate  HIP  solutions  is  to  compare  them  with  the  opti- 
mal mixed-integer  solutions  obtained  using  the  branch-and-boxmd  ap- 
proach. Doing  this  on  a  large  scale,  however,  proved  to  be  cost-prohibitive. 
This  comparison  was  made  with  Small  Twin  Rocks,  a  land  management 
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and  transportation  planning  problem  adapted  from  a  portion  of  the  Twin 
Rocks  area  (Jones  and  others  1986).  The  size  statistics  for  this  model  are 
simimarized  in  table  1.  The  polygons  and  road  network  are  presented  in 
figures  3  and  4,  respectively. 


Table  1 — Size  statistics  for  the  five  models  included  in  this  report 

Small 

Cods  land 

Moose 

Lowinan 

TUJO-RECOYLE 

Catsaorv 

Twin  Rocks 

Creek 

Creek 

area 

8r08 

Polygons 

28 

282 

301 

324 

569 

Road  Links: 

Proposed 

22 

84 

194 

126 

167 

Existing 

2 

64 

79 

104 

67 

Total 

24 

148 

273 

231 

234 

Demand  nodes 

2 

2 

4 

4 

1 

Periods 

3 

3 

2 

3 

3 

Matrix  data: 

Rows 

365 

2.281 

2.150 

3,326 

3,033 

Columns 

260 

1,651 

1,923 

3,515 

4.046 

Integers 

67 

318 

546 

693 

702 

Density^ 

3.1% 

0.8% 

0.4% 

0.4% 

0.4% 

Percentage  of  matrix  coefficients  that  are  nonzero. 


Figure  3— Polygons  for  the  Small  Twin  Rocks 
planning  model. 
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EXISTING  ROADS 


PROPOSED  ROADS 


I     ^/i  MILE 


Figure  4 — Road  network  for  the  Small 
Twin  Rocks  planning  model. 


Six  different  formulations  of  this  model  were  solved  with  both  the  HIP 
process  and  the  branch-and-bound  algorithm  in  the  IBM  Mathematical 
Programming  System  (MPSX).  These  formulations  represent  a  combina- 
tion of  constraints  and  objectives  typical  of  an  actual  planning  problem. 

Table  2  summarizes  the  six  formulations  and  the  resulting  objective 
function  values  for:  (1)  the  first  continuous  LP  solution  in  the  HIP  process, 
(2)  the  HIP  solution,  and  (3)  the  optimal  mixed-integer  solution.  The  objec- 
tive in  formulations  1  through  4  was  to  maximize  present  net  value  (PNV). 
The  side  constraints  in  these  formulations  varied  fi'om  none  (formulation  1), 
to  upper  bounds  of  450  tons  of  additional  sediment  production  in  each 
period  and  lower  boimds  of  3,000  M  bd  ft  of  timber  harvest  in  each  period 
(formulation  4).  The  objective  in  formulation  5  was  to  minimize  the  pro- 
duction of  additional  sediment  while  harvesting  at  least  3,000  M  bd  ft  of 
timber  in  each  period.  The  objective  in  formulation  6  was  to  minimize 
discounted  total  cost  while  meeting  the  same  timber  harvest  targets  as 
formulation  5. 

The  objective  values  of  the  HIP  solutions  were  within  10  percent  of  the 
value  of  the  optimal  mixed-integer  solution  in  all  but  one  case,  formulation 
3.  Thus,  although  the  HIP  process  did  not  find  optimal  mixed-integer 
solutions,  it  did  find  solutions  whose  objective  function  values  closely  ap- 
proximated the  values  for  the  optimal  mixed-integer  solutions. 
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Comparison  With  Plans  developed  via  current  methods  for  tactical  planning  provide,  per- 

Other  Planning  haps,  a  more  useful  basis  for  comparison.  Unlike  comparisons  with  optimal 

Methods  mixed-integer  solutions,  these  comparisons  can  be  done  with  full-scale 

areas.  They  also  measure  the  potential  gains  (or  losses)  in  efficiency  that 
can  be  expected  with  the  HIP  solution  procedure. 

A  common  method  for  analyzing  when  and  where  to  harvest  timber  and 
construct  roads  in  the  Northern  Rockies  utiHzes  professional  judgment  to 
make  the  harvesting  selections  and  a  road  network  cost  minimization 
model,  such  as  NETWORK  (Sessions  n.d.),  to  make  the  roading  selections. 
Plans  developed  via  the  current  HIP  algorithm  were  compared  against  this 
method  on  two  planning  areas  in  the  Northern  Rockies.  The  Copeland 
Creek  area,  on  the  Kootenai  National  Forest  in  Montana,  was  one  of  the 
areas  included  in  an  earher  study  (Jones  and  others  1986).  For  this  area, 
the  HIP  process  developed  a  new  management  alternative,  which  was 
compared  to  a  management  alternative  developed  in  that  earher  study 
using  the  conventional  approach.  The  second  area,  the  Moose  Creek  area, 
is  on  the  Helena  National  Forest.  The  model  for  this  area  was  developed, 
solved,  and  the  results  compared  by  National  Forest  System  personnel 
(Bower  1990).  For  both  areas,  identical  data  and  management  objectives 
were  used  in  the  conventional  and  HIP  approach. 

The  results  of  these  comparisons  are  summarized  in  table  3.  The  objec- 
tive for  the  Copeland  Creek  area  was  to  find  the  combination  of  harvesting 
and  road  building  that  maximizes  present  net  value,  while  maintaining  the 
additional  sedimentation  in  each  of  four  watersheds  below  management- 
specified  upper  bounds.  The  objective  for  the  Moose  Creek  area  was  to  find 
the  combination  of  harvesting  and  road  building  that  maximizes  present 
net  value,  without  exceeding  upper  bounds  on  the  number  of  acres  having 
regeneration  harvests. 


Table  3 — Comparing  HIP  solutions  with  plans  developed  via  a  process  in  which  managers 
choose  the  arrangement  and  timing  of  land  management  activities  and  then  use 
a  network  cost  minimization  model  to  find  the  least-cost  access  routes 


Category 


Copeland  Creek 
area 


Moose  Creek 
area 


Objective 
Side  constraints 


Objective  function 
values: 

(1)  First  continuous  LP 

(2)  HIP 

(3)  Manager  choice 
with  road  cost 
minimization 

Difference  between 
(3)  and  (2)  as  a 
percentage  of  (3) 


Maximize  PNV^ 

Additional  sediment 
in  each  of  four  drainages 
must  not  exceed  specified 
amounts  in  each  period 


4,491 
4.296 

3.014 

43 


Maximize  PNV 

The  number  of  acres 
having  regeneration 
harvests  must  not 
exceed  the  specified 
amounts  in  each  period 


1.665 
1.332 


978 


35 


^Present  net  value  (thousands  of  dollars). 
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In  both  cases  the  HIP  solution  process  developed  plans  that  met  the  con- 
straints and  with  objective  function  values  that  substantially  exceeded  the 
professional  judgment/network  cost  minimization  approach.  For  Copeland 
Creek,  the  increase  was  43  percent;  for  Moose  Creek,  35  percent.  For  the 
Copeland  area,  the  objective  value  for  the  HIP  plan  was  within  5  percent  of 
the  value  for  the  first  continuous  LP,  indicating  the  objective  value  of  this 
solution  closely  approximates  the  optimal  mixed-integer  value.  For  Moose 
Creek  the  objective  value  was  about  21  percent  less  than  the  value  for  the 
first  continuous  LP  solution.  These  results  are  of  approximately  the  same 
magnitude  as  comparisons  made  in  an  earlier  study  (Jones  and  others  1986). 

Comparison  With  The  first  continuous  LP  solution  made  in  the  HIP  process  provides  another 
the  First  basis  for  comparison.  This  LP  solution  is  the  optimal  solution  to  the  con- 

Continuous  LP  tinuous  version  of  a  model  (road  projects  can  assume  any  firactional  value). 

Solution  Restricting  road  projects  to  integer  values  places  an  additional  hmitation 

on  a  model.  This  additional  limitation  will  decrease  the  objective  function 
value  for  maximization  problems,  increase  it  for  minimization  problems,  or 
possibly  have  no  effect  (all  road  projects  have  integer  values  in  the  first 
continuous  LP  solution).  An  additional  limitation  can  never  improve  the 
objective  function  value  of  a  continuous  LP  problem.  The  objective  function 
value  of  the  first  continuous  LP  is,  therefore,  an  upper  bound  for  the  objec- 
tive function  value  of  the  optimal  mixed-integer  solution. 

This  upper  bound  provides  a  convenient,  but  imperfect  yardstick.  It  is 
imperfect  because  it  is  unknown  how  tight  that  bound  is  for  any  one  problem 
(unless,  of  course,  the  optimal  mixed-integer  solution  is  known).  Useful  in- 
formation is  obtained  by  this  comparison  only  when  the  difference  between 
the  HIP  solution  and  the  first  LP  solution  is  small.  A  large  difference  does 
not  necessarily  indicate  a  poor  HIP  solution,  because  it  could  be  reflecting  a 
large  difference  between  the  optimal  MIP  solution  and  the  first  LP  solution. 
For  example,  in  the  Smaill  Twin  Rocks  model  (table  2)  the  difference  between 
the  HIP  and  first  LP  solution  was  quite  large,  while  the  difference  between 
the  HIP  and  optimal  MIP  solution  was  small. 

Table  4  compares  the  objective  function  values  of  15  HIP  solutions  with 
the  respective  values  from  the  first  continuous  LP  solutions.  These  solu- 
tions represent  variations  in  the  constraints,  level  of  constraints,  and  objec- 
tive functions  for  three  full-scale  models.  One  was  the  Moose  Creek  model 
described  earlier.  The  second  model  was  developed  to  analyze  salvage 
harvesting  options  for  the  Lowman  Fire  Complex  on  the  Boise  National 
Forest  in  Idaho  (Bower  1990).  The  third  model  was  developed  for  an  area  on 
the  Lolo  National  Forest  in  Montana  called  TUJO-RECOYLE  (Bower  1990). 
The  size  statistics  for  these  models  are  presented  in  table  1. 

The  HIP  objective  function  values  for  the  Lowman  model  were  within 
5  percent  of  the  value  of  the  first  LP  solution  in  all  but  one  case,  and  within 
10  percent  in  all  cases.  For  these  cases,  the  objective  function  values  of  the 
HIP  solution  closely  approximated  the  mixed-integer  optimums.  For  the 
Moose  Creek  and  TUJO-RECOYLE  models,  the  differences  between  the 
objective  function  value  of  the  first  continuous  LP  solution  and  the  HIP 
solution  were  generally  somewhat  larger.  For  Moose  Creek  the  deviations 
varied  from  24  percent  to  29  percent,  while  for  the  TUJO-RECOYLE  model 
they  varied  from  4  percent  to  21  percent.  These  larger  deviations,  however, 
do  not  necessarily  indicate  poor  solutions  for  the  reasons  discussed  earUer. 
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Table  4 — Comparing  the  objective  values  of  HIP  solutions  made  for  full-scale  models  with  the 

objective  values  of  the  first  continuous  LP  solutions  made  in  the  HIP  process.  (Models 
vary  by  the  objective  function,  and  by  the  level  and  type  of  side  constraints) 


ouiuiion 

Objective  values 

no  la  live 

nuiTiuer 

Mnrlol 

Objective 

First  LP 

HIP 

□eviaiion 

-  -  Thousand  dollars  •  - 

Percent 

1 

Lowman  area 

Max. 

PNV2 

1,686 

1.529 

9.3 

2 

Lowman  area 

Max. 

PNV 

6,381 

6,234 

1.5 

3 

Lowman  area 

Max. 

PNV 

5,788 

5,668 

2.2 

4 

Lowman  area 

Max. 

PNV 

4,602 

4,512 

1.9 

5 

Lowman  area 

Max. 

PNV 

5,308 

5,188 

2.3 

6 

Lowman  area 

Max. 

PNV 

4,896 

4,777 

2.4 

7 

Moose  Creek 

Max. 

PNV 

2,056 

1,504 

26.8 

8 

Moose  Creek 

,  Max. 

PNV 

1,665 

1,227 

26.3 

9 

Moose  Creek 

Max. 

PNV 

1,579 

1.198 

24.1 

10 

Moose  Creek 

Max. 

PNV 

1,932 

1.366 

29.3 

11 

TUJO-RECOYLE 

Max. 

PNV 

-3,802  ■ 

--4,219 

11.0 

12 

TUJO-RECOYLE 

Max. 

VQ|3 

2,616 

2,765 

5.7 

13 

TUJO-RECOYLE 

Max. 

PNV 

-2,038 

-2,276 

11.7 

14 

TUJO-RECOYLE 

Max. 

VOL" 

110 

105 

4.1 

15 

TUJO-RECOYLE 

Max. 

VOL 

110 

86 

21.5 

^The  difference  in  objective  function  values  divided  by  the  value  for  the  first  LP,  multiplied  by  100  to 
express  the  relative  deviation  as  a  percent. 
2PNV  is  present  net  value. 

^VQI  is  a  visual  quality/  index  applying  over  all  periods. 
''VOL  is  timber  harvest  volume  over  all  periods. 


Computer  The  15  full-scale  model  solutions  presented  in  table  4  provide  a  good  basis 

Requirements  for  estimating  the  computer  requirements  of  the  HIP  solution  process.  The 

computation  times  and  costs  associated  with  these  solutions  are  summa- 
rized in  table  5.  Perhaps  the  most  relevant  categories  are  the  CPU  and  I/O 
times  because  they  have  the  most  bearing  in  predicting  performance  on 
other  computer  systems.  (Resource  time  is  actually  a  sum  of  the  other  time 
components  multiplied  by  the  amount  of  memory  used,  and  CC/ER  is  also 
based  on  CPU  time.)  CPU  time  averaged  21.9  minutes  and  I/O  time  aver- 
aged 47.9  minutes.  The  average  total  cost  of  these  solutions,  based  on  the 


Table  5 — Computer  times  and  costs  for  the  15  full-scale  model  HIP  solutions  reported  in  table  3 
run  on  System  C  at  NCC-FC  (UNISYS  Model  1 193) 


 NCC-FC  cost  categories   Total 

Resource  time^     CPUtime^        l/Otlme^       CC/ER  time*  costs 


 Minutes  —  -  

Average                     242.0               21.9               47.9                10.7  $120.07 

Highest                     349.2               37.0               74.1                12.6  198.67 

Lowest                      163.6                13.5               22.5                 8.4  57.37 


^Resource  time,  actually  an  index  for  assessing  a  memory  charge,  is  the  sum  of  CPU,  I/O,  and  CC/ER 
time  multiplied  by  memory  usage. 

^CPU  time  is  central  processing  unit  time. 

^1/0  time  is  time  devoted  to  input-output  devices. 

''CC/ER  is  a  table  look-up  value  based  on  approximate  CPU  time  for  a  typical  user  request. 
^Total  cost  represents  the  composite  measured  in  dollars  using  the  July  1990  price  structure  for  t-priority 
at  NCC-FC. 
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July  1990  T-priority  rates  on  System  C  at  NCC-FC,  was  $120  and  ranged 
from  $57  to  $199.  While  not  small,  these  costs  are  modest  when  compared 
to  the  cost  of  solving  these  models  with  the  mixed-integer  branch-and-bound 
algorithm  available  in  the  Functional  Mathematical  Programming  System 
(Sperry  Corporation  1984)  available  at  NCC-FC.  Earher  attempts  at  solving 
models  of  a  similar  size  using  this  branch-and-boimd  algorithm  typically 
cost  in  excess  of  $500  and  still  did  not  find  the  optimal,  or  in  some  cases 
even  a  feasible,  mixed-integer  solution. 

Below,  the  cost  of  a  HIP  solution  is  expressed  as  a  multiple  of  the  cost  of 
the  first  continuous  LP  solution  made  in  the  HIP  process.  This  solution  is 
always  made  from  the  slack  (beginning)  basis.  Estimates  of  the  resources 
required  to  run  the  HIP  procedure  on  different  computers  can  be  made  by 
multiplying  these  factors  by  the  time  required  to  make  a  continuous  LP 
solution  of  an  IRPM-type  formulation  on  that  computer.  The  multipliers  are: 

HIP  procedure  as  a  multiple 
Category  of  the  first  LP  solution 

Resoiux:e  time  13.7 

CPU  time  8.5 

I/O  time  16.2 

CC/ER  time  24.2 

Total  cost  11.6 

In  the  current  configuration,  the  LP  solver  and  the  program  containing 
the  heuristic  decision  rules  are  written  as  separate  programs.  They  are 
Unked  together  in  a  control  language  loop.  This  means  that  each  time  these 
processes  are  accessed  (t)rpically  seven  to  12  times  for  a  HIP  solution)  the 
internal  arrays  and  variables  must  be  reinitialized  fi*om  data  stored  in  files. 
In  addition,  the  two  processes  communicate  with  each  other  via  files.  There 
is  a  potential  for  substantisd  savings  in  I/O-related  costs  by  integrating  the 
decision  rule  program  and  an  LP  solver  into  one  program. 


CONCLUSIONS 

The  HIP  solution  procedure  was  designed  to  provide  efficient,  feasible 
mixed-integer  solutions  for  the  IRPM-type  formulation  at  a  reasonable  cost. 
To  this  point  we  have  had  good  success  in  achieving  feasible  mixed-integer 
solutions  for  both  test  models  and  full-scale  planning  models.  Costs  for  fiill- 
scale  models  can  be  expected  to  fall  in  the  range  of  $50-$200  in  T-priority, 
given  July  1990  rates  at  NCC-FC.  These  costs  are  modest  compared  to  the 
costs  associated  with  achieving  mixed-integer  solutions  to  these  models 
using  branch-and-bound  algorithms. 

Comparisons  made  with  optimal  mixed-integer  solutions  and  other  plan- 
ning methods  have  been  encouraging.  Plans  developed  using  the  HIP  solu- 
tion procedure  have  been  foimd  to  increase  present  net  value  objectives  by 
as  much  as  40  percent  over  plans  developed  via  a  planning  approach  fre- 
quently used  in  the  Northern  Rockies.  In  comparisons  made  with  optimal 
mixed-integer  solutions,  the  objective  function  value  for  the  HIP  solutions 
was  within  10  percent  of  the  optimal  value  in  five  of  six  cases. 

We  believe  the  HIP  procedure  has  much  potential  for  arranging  and  sched- 
uling land  management  activities  and  road  construction  projects  in  ways 
that  efficiently  meet  specified  management  objectives.  This  efficiency 
should  be  valuable  in  helping  land  managers  meet  the  many  demands  that 
are  placed  on  our  forest  lands. 
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APPENDIX  A:  CAPACITY  ADJUSTMENT  NO.  1 

This  capacity  adjustment  is  done  in  iteration  1.  When  a  problem  is  formu- 
lated, it  is  tjrpical  for  the  capacity  of  the  roads  to  greatly  exceed  the  observed 
quantity  of  traffic  on  those  roads.  As  explained  in  the  main  text  under 
iteration  1,  this  results  in  small  fractional  values  being  assigned  to  the  road 
projects  in  the  first  continuous  LP  solution.  The  objective  of  this  capacity 
adjustment  is  to  perform  a  general  lowering  of  the  capacity  coefficients  so 
that  they  more  closely  match  the  quantities  of  traffic  observed  on  the  road 
network  in  the  first  LP  solution. 

Capacity  adjustment  No.  1  is  made  as  follows: 

where  MC  *    is  the  new  capacity  coefficient  to  be  incorporated  into  the  LP 
matrix,  C*ft  <  is  the  current  capacity  coefficient  for  link  ab  in  period  t,  Cabjt 
is  the  current  capacity  coefficient  for  the  road  standard  immediately  below  k 
(ifk  is  the  lowest  road  standard,  Cabjt  equals  0),  and  AFJ*  is  an  adjustment 
factor  which  is  computed  for  each  road  standard. 

The  first  step  in  computing  AF2*  is  to  calculate  preliminary  adjustment 
factors,  AFlpreiirm  for  each  road  standard: 

AFlprelim  =  the  maximimi  iFab.iJ^ab.p) 
where 

^ab.p  =    the  largest  total  amount  of  traffic  (traffic  is  converted  to  the 
units  of  measure  defined  for  capacity  and  summed  over  all 
traffic  types)  observed  on  link  ab  in  any  one  period.  Super- 
script s  identifies  the  lowest  road  standard  that  can  handle 
that  volume  of  traffic,  and  p  is  the  period  in  which  that  traffic 
was  observed. 

^ab,p  =    the  capacity  coefficient  for  the  road  standard,  s,  and  period, 
p,  associated  with  Fab,p  on  link  ab. 

Separate  adjustment  factors,  AFl^,  are  then  computed  for  each  k  road 
standard.  For  the  highest-capacity  road  standard: 

AFi*  =  AFl^relim 

providing  there  was  at  least  some  observed  traffic  on  this  road  standard.  If 
there  was  no  traffic  on  this  road  standard  (all  traffic  on  the  network  can  be 
carried  by  lower  standards),  AFi*  is  fixed  to  0.15. 

For  road  standards  below  the  highest-capacity  road  standard,  the  adjust- 
ment factor,  AFi*,  is  computed  according  to  the  following  rules.  These  rules 
are  based  on  whether  the  largest  flow,  F^  p,  matches  the  road  standard 
being  processed,  standard  k.  Fabp  is  said  to  match  road  standard  k  if  the 
capacity  of  standard  k  exceeds  Fab,p>  while  the  capacity  for  the  next  lower 
road  standard  (A;  -  1)  is  less  than  Flb,p- 

1.  If  the  largest  flow,  Flb,py  matches  the  road  standard  kik  =  s)  on 
at  least  5  percent  of  the  network  links,  then: 

AF2*  =  the  larger  of  AFi  *^e/tm  or  0.3 

(The  minimum  of  0.3  prevents  excessive  modification.) 
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2.  If  the  largest  flow,  F'ab^p,  matches  road  standard  k(k  =  s)  on  less 
than  5  percent  (or  none)  of  the  network  links,  and  95  percent  or 
more  of  the  flow  on  the  network  is  on  lower  standards,  then: 

AFi*=  0.15 

3.  If  the  largest  flow,  F^  p,  matches  road  standard  k(k  =  s)  on  less 
than  5  percent  (or  none)  of  the  network  links,  but  less  than 

95  percent  of  the  network  flow  is  on  lower  standards,  then: 

AFl''  =  1 

(In  this  case  no  capacity  adjustment  is  made  because  there  is 
insuflficient  information  for  modifying  capacity.) 


APPENDIX  B:  CAPACITY  ADJUSTMENT  NO.  2 

This  capacity  adjustment  is  made  in  iteration  2,  and  in  later  iterations 
(beginning  with  iteration  4)  when  there  is  a  need  for  an  additional  capacity 
adjustment.  The  objective  of  this  adjustment  is  to  further  reduce  the  capac- 
ity coefficients  associated  with  fractional  links  to  promote  larger  fractional 
vsdues  in  subsequent  iterations.  But  care  must  be  taken  to  avoid  reducing 
capacity  so  far  that  traffic  is  redirected  to  inefficient  routes  in.  the  later 
iterations.  This  is  accomplished  by  adjusting  capacity  for  the  road  standard 
matching  the  quantity  of  flow  in  the  cxirrent  LP  solution  to  the  larger  of:  (1) 
the  average  traffic  flow  on  the  network  links  or  (2)  the  observed  flow  on  the 
hnk  being  processed. 

Capacity  adjustment  No.  2  is  based  on  AF2,  a  factor  that  applies  across 
all  road  standards  in  a  model.  The  first  step  in  computing  AF2  is  to  com- 
pute Fib.p  for  each  link: 

Fab,p  =         largest  total  amount  of  traffic  (traffic  is  converted  to  the 
imits  of  measure  defined  for  capacity  and  summed  over  all 
traffic  types)  observed  on  Unk  ab  in  any  one  period. 
Superscript  s  identifies  the  lowest  road  standard  that  can 
handle  that  quantity  of  traffic,  and  p  is  the  period  in  which 
that  traffic  was  observed. 

The  adjustment  factor,        is  computed  as: 

AF2  =    the  larger  of  {0.3,  or  [I  (F^  p/C^  p)VN] 

ab      '  ' 

where  C^,p  is  the  capacity  coefficient  for  the  road  standard  matching  the 
largest  flow  on  link  ah,  which  occurs  in  period p,  and  N  is  the  number  of 
Unks  for  which  a  value  of  F^,p  was  computed  (all  links  having  at  least 
some  traffic  are  included). 

The  rules  for  applying  the  capacity  adjustment  factor  AF2  to  link  ab  vary 
depending  on  whether  the  road  standard  to  be  adjusted  equals,  is  lower,  or 
is  higher  than  the  road  standard  matching  the  highest  traffic  volume, 
Fib.pt  on  the  link.  The  matching  road  standard  is  the  lowest  road  standard 
that  will  handle  the  highest  traffic  volume  on  the  Unk. 


26 


1.  If  the  road  standard  being  adjusted  is  standard  s,  the  matching 
standard  for  the  link,  capacity  for  each  period  (t)  is  adjusted 
according  to  the  following  rules: 

(i)  If  (AF2XCofe,p)  ^  Flb.p,  then  the  capacity  for  standard  s  is  ad- 
justed for  each  period  as  follows: 

MC^b.t  =  (AF2XC^b,t)  (for  each  t) 

(ii)  If  iAF2XClb,p)  <  Fab.p,  then  the  adjustment  for  each  period  is: 

MC^b.t  =  (F^b,pXC^b.t  IC^b,p)  (for  each  t) 

MCab.t  is  the  new  capacity  coefficient  to  be  incorporated  into  the 
LP  matrix  for  standard  s  in  period    Cib.t  is  the  previous 
capacity  coefficient  for  standard  s  in  period  t{t  =  1,  2,  etc.),  and 
Cab,p  is  the  previous  capacity  coefficient  for  standard  s  in 
period  p,  the  period  of  highest  observed  flow  on  Unk  ab.  The  ratio 
(Cab,t/Cab,p)  is  included  to  adjust  F^b.p  when  the  capacity  for  road 
standard  s  varies  across  planning  periods.  This  occurs  when  the 
periods  are  of  unequal  length. 

2.  If  the  road  standard  being  adjusted  is  higher  than  standard  s  (the 
lowest  standard  that  can  handle  the  observed  flow  on  Unk  ab) 
capacity  is  adjusted  as  follows: 

k        k-i  k  k-i 

MCab.l  =  Cab.i  +  (AF2)(Cab,t  -  Cab,t) 

where  Cab.t  is  the  previous  capacity  coefficient  for  standard  k  on 
hnk  ab  in  period  t,  and  Cab}  is  the  previous  capacity  coefficient  in 
period  t  for  the  next  lower  road  standard,  k-1. 

3.  If  the  road  standard  being  adjusted  is  lower  than  standard  s  (the 
lowest  standard  that  can  handle  the  observed  flow  on  Unk  ab) 
capacity  is  not  adjusted.  There  is  no  reason  to  adjust  capacity 
downward  because  the  flow  on  hnk  ab  exceeds  the  capacity  for  this 
road  standard. 

4.  If  Unk  ab  carries  no  traffic  in  the  current  LP  solution,  the 
capacity  coefficients  are  modified  as  follows  for  all  the  standards 
defined  for  that  link: 

MCab,t  =  Cab,i+  (AF2XCab.t  "  Cab.l) 

where  all  terms  are  as  previously  defined.  Ifk  is  the  lowest  road 
standard  defined  for  Unk  ab,  then  Cab}  equals  0. 


APPENDIX  C:  ROUNDING  PROCEDURE 

The  rounding  procedure  consists  of  three  major  steps:  (1)  develop  a  list  of 
candidate  road  projects  for  rounding,  (2)  compute  an  index  for  ranking  the 
candidate  projects,  (3)  choose  the  projects  to  be  rounded  to  1  from  the  set  of 
candidate  projects  connecting  to  either  a  final  demand  node  or  to  road 
projects  that  were  previously  roimded  to  1. 
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Develop  a  List  of 
Candidate  Road 
Projects  for 
Rounding 


One  candidate  project  is  selected  for  each  fractional  link  that  carried 
traffic  in  the  current  LP  solution,  or  was  associated  with  a  resource  project 
that  contributed  flow  to  the  road  network.  (Association  with  resource 
projects  is  through  resource  project-to-road  triggers.  See  equation  8  in  the 
general  formulation.)  First,  a  preferred  period  and  road  standard  are 
selected.  Then,  a  candidate  project  is  selected  after  comparing  the  preferred 
period  and  standard  selections  with  the  feasible  projects  available  for  that 
Hnk. 

Preferred  Period — The  preferred  period  for  link  ab,  Pab,  is  the  earliest 
period  in  which  there  is  a  significant  quantity  of  traffic  on  link  ab.  Pab  is 
the  earliest  period  in  which: 

Tab,t  ^  MFLOWt 

where  Tab,t  is  the  total  traffic  on  link  ab  in  period  t  (each  traffic  tj^je  is 
converted  to  capacity  units  and  simimed)  and  MFLOWt  equals  the  smaller 
nonzero  value  of:  MFLOWlt  or  MFL0W2t.  MFLOWlt,  a  minimum  flow 
criterion  based  on  the  average  network  flow,  is  calculated  as: 

MFLOWlt  =  (0.05)  (Z  Tab.tVN 

where  Tab,t  is  as  defined  above,  and  N  is  the  number  of  links  having  non- 
zero flow  in  period  t.  MFL0W2t  is  a  minimum  flow  criterion  based  on  the 
flow  loaded  onto  the  network  by  resource  projects.  For  all  periods  except 
the  last  period: 

MFLOW2t  =  {V2XFIXRPt) 

where  FLKRPt  is  the  smallest  quantity  of  traffic  loaded  onto  the  network  in 
period  t  from  any  one  resource  project  that  was  preselected  (set  into  solu- 
tion) by  the  user.  If  a  resource  project  creates  flow  in  more  than  one  period, 
however,  only  the  initial-period  traffic  is  considered  in  calculating  MFL0W2t . 
For  the  last  period: 

MFL0W2t  =  {V2tRPt) 

where  RPt  is  the  smallest  quantity  of  traffic  loaded  onto  the  network  in  the 
last  period  by  any  resoiu'ce  project  (whether  fixed  by  the  user  or  not)  whose 
first  traffic  loading  occurs  in  the  last  period. 

MFL0W2t  is  designed  to  ensure  that  period  selection  is  consistent  with 
providing  access  to  necessary  resource  projects.  In  all  but  the  last  period, 
necessary  projects  are  those  preselected  by  the  user.  For  the  last  period,  all 
resource  projects  are  considered  necessary  to  accommodate  the  tranship- 
ment formulation  where  at  least  one  harvesting  option  must  be  selected 
from  each  polygon  (in  the  transhipment  formulation  equation  4  is  modified 
from  an  upper-bound  constraint  to  an  equahty). 

Preferred  Road  Standard — The  preferred  road  standard  for  link  ab, 
Sab,  is  the  lowest  road  standard  whose  capacity  exceeds  both  (1)  the  largest 
total  traffic  flow  observed  on  that  link  in  any  one  period  and  (2)  the  largest 
total  traffic  flow  loaded  onto  the  network  in  any  one  period  from  the  re- 
source projects  associated  with  link  ab  (through  resource  project-to-road 
triggers).  These  total  flows  include  all  traffic  types  converted  to  capacity 
units.  If  the  observed  quantity  of  traffic  in  any  one  period  exceeds  the 
capacity  for  the  highest  capacity  road  standard  defined  for  that  link,  the 
highest  capacity  standard  is  selected.  (After  the  road  project  is  roimded 
to  1,  capacity  becomes  an  upper  bound  on  the  total  quantity  of  traffic.) 
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Candidate  Selection — The  selection  of  the  candidate  project  for  link  ab 
is  based  on  the  following  rules: 

1.  If  the  preferred  road  standard,  Sofe,  represents  an  existing  road  requir- 
ing no  reconstruction  (there  is  no  cost  associated  with  the  road  project), 
then  assign  a  value  of  1  to  the  road  project  on  link  ab  having  that  standard 
in  the  earliest  period  available.  This  road  already  exists  on  the  ground,  and 
the  decision  rules  indicate  that  it  should  not  be  replaced. 

2.  If  the  preferred  road  standard.  Sab,  is  not  an  existing  road,  but  is 
available  as  a  proposed  project  in  the  preferred  period,  Pab,  then  that 
project  is  selected  as  the  candidate  for  link  ab,  providing  that  selection  is 
not  precluded  by  a  project  preselected  by  a  user.  Preselected  projects  in- 
clude road  construction  projects  for  Unk  ab  and  resource  projects  that 
trigger  other  road  projects  for  Unk  ab  that  were  incorporated  into  the  solu- 
tion by  the  user. 

3.  If  the  preferred  road  standard.  Sab,  is  not  an  existing  road  and  is  not 
available  as  a  proposed  project  in  the  preferred  period,  Pab  (because  that 
project  is  infeasible  as  a  result  of  a  road  or  resource  project  preselected  by 
the  user,  or  because  a  project  of  standard  Sab  was  not  defined  for  period 
Pab),  then  an  alternate  project  is  selected.  The  first  priority  is  to  find  a 
project  with  a  different  road  standard  that  has  the  same  implementation 
period.  If  the  quantity  of  traffic  on  that  link  approaches  the  capacity  of  the 
selected  road  standard,  projects  with  higher  road  standards  are  preferred. 
Otherwise  projects  with  lower  road  standards  are  preferred.  If  no  accept- 
able project  during  the  same  implementation  period  is  found,  projects  with 
earlier  implementation  periods  are  investigated.  There  must  be  at  least  one 
project  feasible  for  the  observed  traffic  in  the  current  solution  to  have 
occurred.  The  details  of  this  selection  process  follow.  Let  Fab,p  represent 
the  largest  traffic  flow  for  link  ab;  the  subscript  p  indicates  the  period  in 
which  the  largest  flow  occiured.  Let  Cab.p  represent  the  capacity  for  the 
preferred  road  standard  in  period  p,  and  Cab,p  represent  the  capacity  in 
period  p  for  the  next  lower  road  standard. 

If'     ^ah.p  ^  Cab.p  +  0  .75iCab,p  -  Cab.p) 

then  select  the  project  with  the  next  higher  road  standard,  providing  it  is 
feasible  with  regard  to  any  road  or  resource  projects  associated  with  link  ab 
that  were  preselected  by  the  user.  If  no  feasible  higher  standards  are  avail- 
able, select  the  project  with  next  lower  standard,  providing  it  is  feasible.  If 
no  feasible  road  standards  are  available  in  the  preferred  period,  repeat  this 
process  for  period  Pafc-1,  then  Pafe-2,  and  so  on.  There  must  be  an  available 
road  project  in  one  of  these  periods  for  traffic  to  have  appeared  on  link  ab  in 
time  Pab- 

If'     I^ab.p  <  Cab.p  +  0  •'^^iCab.p  -  Cab.p) 

then  select  the  project  with  the  next  lower  road  standard,  providing  it  is 
feasible  with  regard  to  any  road  or  resource  projects  associated  with  link  ab 
that  were  preselected  by  the  user.  If  no  feasible  lower  standards  are  avail- 
able, select  the  project  with  next  higher  standard,  providing  it  is  feasible.  If 
no  feasible  road  standards  are  available  in  the  preferred  period,  repeat  this 
process  for  period  Pab-^,  then  Po6-2  and  so  on.  There  must  be  an  available 
road  project  in  one  of  these  periods  for  traffic  to  have  appeared  on  link  ab  in 
time  Pab- 
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Calculate  Index  The  next  step  is  to  calculate  an  index  for  ranking  the  candidate  road 

Values  for  the  projects.  The  smaller  the  value  of  this  index,  the  more  preferred  the  road 

Candidate  project.  The  index  for  road  project  j  {Indexj)  is  calculated  as  follows: 

Projects  Indexj  =  3(W/,)  -  WFj  +  WCj 

where  WIj  measures  the  relative  importance  of  road  project  j  in  the  side 
constraints  (the  larger  the  value  the  less  desirable  the  index),  WFj  measures 
the  relative  amount  of  traffic  flow  on  link  ab  for  which  project  j  is  an  option 
(the  more  the  traffic  the  better  the  index),  and  WC/ measures  the  relative 
cost  of  road  project  j  in  the  objective  function  (the  larger  the  cost  the  less 
desirable  the  index).  The  computations  for  each  component  are  described  in 
detail  below. 

The  WIj  Component — WIj  is  the  relative  importance  of  candidate  project 
y  in  the  side  constraints,  weighted  by  1  minus  the  activity  level  for  road 
project  j  in  the  current  LP  solution.  It  is  the  sum  of  the  relative  impacts 
(/MPy)  of  project  J  in  each  of  the  side  constraints: 

n 

WIj=I.IMPij 
i=l 

For  constraints  having  only  upper  bounds 

IMP^j  =[il-Aj)a.jVH. 
For  constraints  having  only  lower  bounds 

IMPij  =  [a-Aj)(-l)aijVLi 
For  constraints  having  both  upper  and  lower  bounds 

/MPy  =  the  larger  of:  [(1 -A,)ay]/i/i 

or 

[(l-AjX-DayVLi 

For  equality  constraints 

IMPij  =  \a-Aj)aij\  lEi 
where 

=    the  coefficient  for  road  project  j  in  side  constraint  i 
Aj     =    the  activity  level  (solution  value)  for  road  project^  in  the 

current  LP  solution 
Hi     =    the  space  available  in  side  constraint  i,  which  is  a  "less  than" 
constraint.  It  is  the  right-hand-side  value  (upper  bound)  of 
constraint  i,  minus  the  activity  level  for  constraint  i  (including 
the  deviation  variables),  plus  the  upper  bound  for  the  lower- 
cost  deviation  variable  for  constraint  i,  plus  0.001.  The  last 
term  is  added  to  avoid  the  possibility  of  dividing  by  0. 
Li      =    the  space  available  in  side  constraint  i,  which  is  a  "greater 

than"  constraint.  It  is  the  activity  level  for  constraint  i,  minus 
the  right-hand-side  value  (lower  bound)  of  constraint  i,  plus  the 
upper  bound  for  the  lower-cost  deviation  variable  for  constraint 
i,  plus  0.001. 

Ei     =    a  proxy  for  room  remaining  in  side  constraint  i,  which  is  an 
equahty  constraint.  It  is  the  right-hand-side  value,  plus  the 
upper  bound  for  the  lower-cost  deviation  variables  (these 
bounds  are  the  same  number)  multiphed  by  0.15.  This 
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r 

^ab.t  = 


calculation  attempts  to  produce  a  proxy  having  the  same 
magnitude  relative  to  its  right-hand-side  value  as  Lj  and  Hi 
have  to  their  respective  bounds.  Ei  is  needed  because  an 
equahty  constraint  will  always  be  exactly  satisfied  (if  the  LP 
solution  was  feasible)  and  therefore  will  never  have  space 
remaining  as  do  upper  and  lower  bound  constraints. 
n       =    the  number  of  side  constraints. 

The  larger  the  value  for  WIj,  the  less  desirable  is  the  index.  The  contribu- 
tion of  Oy  to  the  index  is  inversely  proportional  to  activity  level  of  project  j  in 
the  current  LP  solution  as  a  result  of  the  1  -Aj  weights.  Thus,  the  impor- 
tant road  projects  in  the  continuous  LP  (as  measured  by  the  magnitude  of 
the  activity  level)  have  less  penalty  for  side-constraint  impact  than  do  road 
projects  having  lower  fi^actional  values. 

The  WFj  Component — This  component  measures  the  amount  of  traffic 
flow  observed  on  link  ab  relative  to  the  average  traffic  flow  on  the  network 
in  the  current  LP  solution.  The  larger  this  component,  the  lower  (more 
desirable)  the  index.  It  is  computed  as: 

T  T 

WF,-  =    [ir„fe.,,M(  I      I  W/M] 
t=\   ^     abeL  t  =  1 

where 

^o6,.<  =    total  amount  of  traffic  flow  (expressed  in  capacity  units)  on 
Hnk  abj  (the  Unk  associated  with  road  project  f)  in  period  t 
total  amount  of  traffic  flow  (expressed  in  capacity  units)  on 
hnk  ab  in  period  t 

M     =    total  number  of  nonzero  values  for  Tab,t  across  all  links  on  the 

network 
T      =    number  of  periods 
L      =    number  of  links. 

The  WCj  Component — This  component  is  the  relative  cost  of  the  road 
project,  weighted  by  1  minus  the  activity  level  for  road  project  j  in  the  cur- 
rent LP  solution.  It  is  computed  as  follows  for  maximization  problems: 

J 

WCj=[a-AjK-l)Cj]/[I,  \Cj\/N] 

J=l 

and  for  minimization  problems: 

J 

WCj  =  [(l-Aj)Cj]/[I.  \Cj\/N] 
7=1 

where 

Cj     =    the  objective  function  coefficient  for  project  j 

Aj      =    the  activity  level  for  road  project  j  in  the  current  LP  solution 

N      =    the  total  number  of  nonzero  Cj  values 

J      =    the  total  number  of  candidate  road  projects. 

The  greater  this  component,  the  larger  (less  desirable)  the  index.  The 
weights  (1  -  Aj)  make  the  contribution  of    to  the  index  inversely  propor- 
tional to  the  activity  level  for  project  j  in  the  current  LP  solution.  Thus,  the 
more  important  projects  (as  measured  by  the  size  of  their  fi-actional  value) 
have  less  penalty  for  construction  cost. 
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The  "tree-growth"  rounding  procedure  derives  its  name  from  the  order  in 
which  road  projects  are  roimded.  Only  road  projects  that  connect  either 
directly  to  a  demand  node,  or  indirectly  through  previously  rounded 
projects  can  be  rounded  to  1.  As  more  road  projects  are  rounded  to  1,  the 
growth  of  the  "developed  portion"  of  the  network  resembles  tree  growth. 
Unlike  most  trees,  however,  the  branches  originating  from  one  demand 
node  can  connect  (grow  together)  with  branches  from  other  demand  nodes. 
This  method  of  selecting  road  projects  ensures  that  the  "developed  portion" 
of  the  network  is  comprised  of  one  or  more  continuous  networks,  each  of 
which  connect  to  at  least  one  final  demand  node.  The  following  steps 
describe  the  rounding  process  that  is  applied  in  each  iteration,  beginning 
with  iteration  3. 

1.  Initialize  the  "developed  portion**  of  the  road  network — ^The 
"developed  portion"  of  the  network  consists  of  road  projects  having  integer 
values  which  are  available  for  carrying  traffic  to  the  final  demand  nodes. 
At  the  beginning  of  an  iteration,  the  "developed  portion"  of  the  network  is 
comprised  of  three  parts.  First,  it  includes  the  road  projects  rounded  to  1 
in  the  previous  iterations  (in  iteration  3,  this  consists  only  of  the  demand 
nodes).  Second,  the  candidate  road  projects  representing  existing  roads  are 
rounded  to  1  and  added  to  the  "developed  portion"  of  the  network.  Third, 
any  preselected  road  projects  (user  set  them  to  1)  are  added  to  the  "devel- 
oped portion"  of  the  network. 

2.  Build  the  list  of  "connecting  projects** — The  "connecting  projects" 
are  the  candidate  projects  that  connect  both  spatially  and  temporally  with 
the  "developed  portion"  of  the  network.  A  candidate  project  connects  spa- 
tially if  it  shares  a  node  with  one  or  more  links  included  in  the  "developed 
portion"  of  the  network.  That  project  connects  temporally  if  its  period  of 
implementation  is  greater  than  or  equal  to  the  period  of  implementation  for 
at  least  one  road  construction  project  in  the  "developed  portion"  of  the 
network  with  which  the  common  node  is  shared.  In  other  words,  the  road 
network  must  be  developed  to  the  common  node  by  the  period  of  implemen- 
tation of  a  candidate  project  for  that  project  to  be  feasible  temporally. 

3.  Check  the  feasibility  of  the  "connecting  project**  with  the  low- 
est index  value — First  identify  the  "connecting  project"  having  the  lowest 
(best)  index  value.  If  there  are  no  remaining  "connecting  projects,"  proceed 
to  step  8.  Let  Xcur  represent  the  "connecting  project"  having  the  lowest 
(best)  index  value.  The  feasibility  of  rounding  Xcur  to  1  is  checked  for  each 
side  constraint.  This  is  accomplished  for  constraint  i  by  comparing  the 
impact  of  Xcur  plus  the  cumulative  impacts  from  previous  rounding  deci- 
sions made  in  the  current  iteration  against  the  "allotted  space"  for  con- 
straint i.  "Allotted  space"  is  the  amount  of  impact  irom  rounding  decisions 
that  is  permitted  in  the  current  iteration.  "Allotted  space"  intentionally 
underestimates  the  estimated  "actual  space"  (the  estimated  actual  amount 
of  space  available  for  absorbing  impacts  from  rounded  road  projects)  to 
avoid  making  too  many  decisions  in  any  one  iteration.  (See  appendix  D  for 
details  regarding  this  feasibility  check.) 

(a)  If  the  cumulative  impact  from  rounding  Xc^r  is  within  the  "allotted 
space"  for  all  side  constraints,  go  to  step  5. 

(b)  If  the  cumulative  impact  from  rounding  Xcur  exceeds  the  "allotted 
space"  for  any  side  constraint,  go  to  step  4. 
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4.  Check  the  feasibility  of  Xcur  against  the  '^actual  space**  available 
in  each  side  constraint — In  this  step  the  impacts  of  Xcur  plus  the  cumula- 
tive impacts  from  the  previous  rounding  decisions  made  in  the  current 
iteration  are  compared  against  the  "actual  space"  available  in  each  side 
constraint. 

(a)  If  the  cumulative  impact  associated  with  rounding  Xc^^  is  less  than  the 
"actual  space"  for  all  side  constraints,  then: 

(i)  if  an  alternate  rounding  choice  (Xaii)  lias  not  yet  been  selected  since 
last  updating  the  hst  of  connecting  projects,  let  Xait  =Xcur, 

(ii)  otherwise,  maintain  the  previously  selected  Xq/^  (Xq^/ remains  a 
candidate  for  rounding,  and  will  eventually  be  rounded  providing  its 
index  value  proves  to  be  substantially  better  than  the  next  best 
feasible  project.)  Retimi  to  step  3  to  process  the  "connecting  project" 
with  the  next  best  index. 

(b)  If  the  cumulative  impact  associated  with  roimding  Xcur  exceeds  the 
"actual  space"  for  any  one  side  constraint,  remove  Xcur  from  the 
"connecting  project"  Ust  and  add  it  to  the  "infeasible  connecting  project" 
Hst,  then  return  to  step  3  to  process  the  next  "connecting  project." 

5.  Make  a  rounding  decision — This  step  is  reached  if  the  impacts  from 
Xcur  were  found  by  step  3  to  be  feasible  for  all  side  constraints  (impacts  are 
within  "allotted  space").  The  index  value  for  Xcur,  (INDEXcur)  is  compared 
with  the  index  value  for  the  alternate  project  (INDEXaii)  to  determine  which 
project  should  be  rounded  to  1.  If  the  difference  in  index  value  is  less  than 
or  equal  to  0.5,  Xcur 

is  selected.  In  this  case  Xait  is  not  sufficiently  better  to 
warrant  exceeding  the  threshold  of  acceptable  impact  for  at  least  one  side 
constraint.  If  the  difference  in  index  value  exceeds  0.5,  Xait  is  selected  over 
Xcur  even  though  it  will  have  a  significant  impact  on  at  least  one  side  con- 
straint. The  details  of  this  step  follow: 

(a)  If:  INDEXcur  -  INDEXait  <  0.5,  then  round  Xcur  to  1.  Next,  go  to  step  6. 

(b)  If:  INDEXcur  -  INDEXait  >  0.5,  then  round  Xait  to  1.  Next,  go  to  step  10 
(no  further  rounding  decisions  will  be  made  in  the  current  iteration 
because  Xait  exceeds  the  "allotted  space"  for  at  least  one  side 
constraint). 

6.  Update  various  lists — After  a  road  project  has  been  rounded  to  1, 
various  lists  must  be  updated.  First,  the  rounded  project  is  added  to  the 
"developed  portion"  of  the  network.  Second,  the  impacts  of  the  rounded 
project  are  added  to  the  cumulative  totals  for  the  side  constraints.  Third, 
the  list  of  "connecting  projects"  is  updated  as  follows: 

(a)  The  rounded  project  is  removed. 

(b)  The  candidate  projects  associated  with  fractional  links  connecting  to  the 
rounded  project  are  added. 

(c)  The  "infeasible  connecting  projects"  and  Xait  are  added.  There  is  a 
possibihty  that  the  rounded  project  created  more  "actual  space"  (rather 
than  less)  in  critical  side  constraints.  So,  the  previously  infeasible 
projects  are  added  back  to  "connecting  project"  list  because  they  may 
now  be  feasible. 

Proceed  to  step  7. 
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7.  Has  the  ipn-gimiim  number  of  rounding  decisions  for  the  cur- 
rent iteration  been  reached? — The  maximuin  number  of  projects  that 
can  be  rounded  in  iterations  3  and  4  equals  50  percent  of  the  original  candi- 
date projects  plus  1.  This  is  increased  to  75  percent  plus  1  beginning  with 
iteration  5.  These  restrictions  represent  a  compromise.  Experience  sug- 
gests that  better  solutions  are  obtained  when  a  small  number  of  rounding 
decisions  are  made  in  each  iteration.  This  permits  the  ensuing  LP  solution 
to  react  to  these  rounding  decisions  prior  to  making  additional  decisions. 
On  the  other  hand,  making  as  many  decisions  as  possible  in  each  iteration 
reduces  the  number  of  iterations,  thereby  reducing  the  amount  of  computer 
resources  required  for  a  HIP  solution. 

(a)  If  the  maximum  number  of  roimding  decisions  has  not  been  reached, 
go  back  to  step  3. 

(b)  If  the  maximxmi  number  of  rounding  decisions  has  been  reached, 
proceed  to  step  10. 

8.  Round  the  alternative  project  iXait),  if  one  has  been  identified 
since  updating  the  list  of  "connecting  projects** — This  step  is  reached 
(from  step  3)  if  none  of  the  remaining  "connecting  projects"  is  within  the 
"allotted  space,"  but  space  for  absorbing  impacts  remains  in  the  side  con- 
straints and  the  maximimi  nimiber  of  roimding  decisions  for  the  current 
iteration  has  not  yet  been  reached. 

(a)  If  there  is  an  alternative  project,  Xait,  round  this  project  to  1.  Then 
proceed  to  step  10. 

(b)  If  there  is  no  Xait  available  (the  impact  of  each  connecting  project  plus 
the  cumulative  impact  of  all  previously  rounded  projects  exceeds  the 
"actual  space"  in  each  of  the  side  constraints),  proceed  to  step  9. 

9.  Have  any  rounding  decisions  been  made  in  the  current  itera- 
tion?— This  includes  modifications  made  in  previous  rounding  decisions 
(road  standard  changes,  etc.)  as  well  as  whether  any  projects  have  been 
roimded  to  1  in  the  current  iteration. 

(a)  If  at  least  one  rounding  decision  has  been  made  in  the  current 
iteration,  proceed  to  step  10. 

(b)  If  no  rounding  decisions  have  been  made  in  the  current  iteration,  then 
the  impacts  from  each  of  the  "connecting  projects"  exceed  the  "actual 
space"  in  each  of  the  side  constraints.  Yet,  fi*actional  projects  remain. 
Something  must  be  done  to  avoid  getting  the  same  result  in  the  next 
iteration.  Implement  the  following  rules  and  then  proceed  to  step  10. 

(i)  If  either  the  objective  function  value  or  the  list  of  candidate  projects 
differ  from  the  previous  iteration,  then  each  of  the  "connecting 
projects"  are  set  to  0.  In  addition,  the  road  projects  of  the  same  road 
standard  applying  to  earUer  periods  on  those  links  are  also  set  to  0. 

(ii)  If  both  the  objective  function  value  and  list  of  candidate  projects  are 
the  same  as  the  previous  iteration,  then  set  to  0  all  of  the  candidate 
projects  (whether  they  connect  or  not),  plus  the  road  projects  of  the 
same  road  standard  applying  to  earHer  periods  on  those  links. 

10.  End  the  rounding  procedure  in  the  current  iteration. 
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APPENDIX  D:  FEASIBILITY  CHECKING 

In  the  HIP  process  the  feasibility  of  rounding  option  j  in  side  constraint  i  is 
checked  against  either  "actual  space"  or  "allotted  space."  These  checks  are 
made  in  the  rounding  procedure  (see  the  discussion  for  iteration  3  in  the 
main  report,  and/or  appendix  C)  and  in  modifications  made  to  previous 
rounding  decisions  (see  the  discussion  for  iterations  4,  5,  and  6  in  the  main 
report). 


Comparing 
Against  the 
"Actual  Space" 


"Actual  space"  measures  the  total  amount  of  impact  that  constraint  i  can 
absorb  from  rounding  decisions  in  the  current  iteration.  The  feasibility  of 
rounding  option  j  with  regard  to  the  "actual  space"  in  side  constraint  i  is 
checked  by  comparing  its  impact  {IMPACT ij)  plus  the  cumulative  impacts  on 
constraint  i  of  previous  rounding  decisions  in  the  current  iteration  (CUMi), 
against  the  "actual  space"  in  constraint  i  (SPACE i). 

IMPACTij  is  calculated  for  the  various  types  of  rounding  options  as  follows: 

1.  The  impact  on  constraint  i  of  rounding  road  project  j  to  1  is: 

IMPACTij  =  aij 
where  a-j  is  the  coefficient  of  project  j  in  side  constraint  i. 

2.  The  impact  on  constraint  i  of  a  change  in  road  standard  or  postponing 
construction  in  the  "developed  portion"  of  the  network  is: 

IMPACTij  =  (ay(new))  -  (OyCold)) 

where  a,^new)         coefficient  in  constraint  i  for  the  road  project  replacing 
the  previously  selected  project,  whose  coefficient  in  constraint  i  is  ay(oid)* 

3.  The  impact  on  constraint  i  of  closing  a  link  (by  setting  the  previously 
rounded  road  project  to  0)  is: 

IMPACTij  =  -(a,;^oid)) 

where  Oj^old)  ^®      coefficient  in  constraint  i  of  the  previously  selected 
project. 

CUMi  is  the  cumulative  total  impact  on  constraint  i  of  the  roimding  deci- 
sions made  in  the  current  iteration.  The  impact  of  rounding  option  j  is 
added  to  CUMi  when  rounding  option  j  is  selected  and  incorporated  in  the 
LP  matrix. 

Feasibility  for  Upper  Bound  Constraints — The  available  space  in  side 
constraint  i  for  absorbing  impacts  from  heuristic  rounding  decisions, 
SPACE i,  is  calculated  at  the  beginning  of  a  heuristic  iteration  (immediately 
following  a  LP  solution).  For  upper  bound  constraints  the  calculation  is: 


SPACEi  =  RHSi  -  [Ai  +  DiL  +  Di^-Z ai/Rj)] 


(D-1) 


where 


RHSi  = 
Ai  = 


- 


the  right-hand-side  value  for  side  constraint  i 

the  activity  level  for  side  constraint  i  (includes  the  effects  of 

the  deviation  variables) 

the  activity  level  (solution  value)  for  the  lower-cost  deviation 
variable  associated  with  the  upper  bound  placed  on  constraint  i 
the  activity  level  (solution  value)  for  the  higher-cost  deviation 
variable  associated  with  the  upper  bound  placed  on  constraint  i 
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Oij     =    the  coefficient  for  road  project  j  in  side  constraint  i 
Rj     =    the  activity  level  (solution  value)  for  road  project  j,  which  is  in 
the  current  LP  solution  at  some  fractional  value. 

For  upper  bound  constraints,  a  positive  value  for  SPACEi  indicates  there  is 
"actual  space"  available  for  absorbing  impacts  from  rounding  decisions,  a 
negative  value  indicates  that  the  upper  bound  on  constraint  i  was  violated 
in  the  current  LP  solution  (there  is,  in  fact,  negative  "actual  space"),  and  a 
value  of  zero  indicates  the  constraint  was  at  its  bound  (no  space  is  available). 

Rounding  options  with  negative  IMPACTij  values  decrease  the  cumulative 
total  impact  from  rounding  decisions  {CUMi).  These  options  are,  therefore, 
always  feasible,  regardless  of  the  value  of  SPACEi  or  CUMi. 

Rounding  options  with  positive  values  for  IMPACTij  increase  the  cumula- 
tive total  impact  from  rounding  decisions  (CUMi).  These  options  are  fea- 
sible if 

IMPACTij  +  CUMi  <  SPACEi 

The  cumulative  total  impact  of  rounding  decisions  on  side  constraint  i, 
CUMi,  must  be  updated  following  each  rounding  decision.  This  is  done  by 
adding  IMPACTij  for  the  selected  rounding  option  to  CUMi. 

Feasibility  for  Lower  Bound  Constraints — The  available  space  with 
regard  to  a  lower  bound  placed  on  side  constraint  i  is  calculated  as  follows: 

SPACEi  =  RHSi  -  [Ai  -DiL-DiH-'^  ai/Rj)]  (D-2) 

where  D/l  is  the  activity  level  for  the  lower-cost  deviation  variable  associ- 
ated with  the  lower  bound  on  constraint  i,  Djjj  is  the  activity  level  for  the 
higher-cost  deviation  variable  associated  with  the  lower  bound  on  con- 
straint i,  and  all  other  terms  are  as  previously  defined. 

A  negative  value  for  SPACEi  indicates  there  is  "actual  space"  available  for 
absorbing  impacts  from  rounding  decisions,  a  positive  value  indicates  that 
the  lower  bound  on  constraint  i  was  violated  in  the  current  LP  solution 
(there  is,  in  fact,  negative  "actual  space"),  and  a  value  of  0  indicates  con- 
straint i  was  at  its  lower  bound  (no  space  is  available). 

Impacts  for  rounding  options  are  calculated  the  same  way  as  for  upper- 
bound  constraints,  but  the  test  for  identifying  feasible  rounding  options  is 
reversed.  Rounding  options  with  positive  /MPACTy  values  increase  CUMi, 
thereby  moving  it  away  from  SPACEi.  These  rounding  options  are,  there- 
fore, always  feasible  for  lower  bound  constraints,  regardless  of  the  values  of 
SPACEi  or  CUMi. 

Rounding  options  with  negative  values  for  IMPACTij  decrease  CUMi 
thereby  moving  it  closer  to  the  "actual  space"  for  constraint  i.  Rounding 
options  with  a  negative  IMPACTij  are  feasible  if 

IMPACTij  +  CUMi  >  SPACEi 

The  cumulative  total  impact  of  rounding  decisions  on  side  constraint  i, 
CUMi,  must  be  updated  following  each  rounding  decision.  This  is  done  by 
adding /MPACTy  for  the  selected  rounding  option  to  CUMi. 

Feasibility  for  Equality  Constraints — Equahty  constraints  have 
upper  and  lower  bounds  set  to  the  same  right-hand-side  value.  Fesisibility 
with  regard  to  both  bounds  is  checked.  A  rounding  decision  is  considered 
infeasible  if  the  calculations  for  either  the  upper  or  lower  bounds  indicate 
an  infeasibihty. 
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Feasibility  for  Constraints  Having  Ranges — Constraints  for  which 
"ranges'*  have  been  specified  contain  both  lower  and  upper  bounds.  Feasibil- 
ity with  regard  to  both  bounds  is  checked.  A  rounding  decision  is  considered 
infeasible  if  the  calculations  for  either  the  upper  or  lower  bounds  indicate  an 
infeasibihty. 

The  "allotted  space"  is  computed  by  adjusting  "actual  space."  When  road 
construction  plays  a  significant  role  in  side  constraint  i  (10  percent  or  more 
of  the  impacts  in  the  constraint  are  associated  with  road  construction), 
"actual  space"  is  adjusted  downward  to  compute  "allotted  space."  In  these 
instances  "allotted  space"  is  a  more  conservative  assessment  of  the  amount 
of  rounding-decision  impacts  to  permit  on  constraint  i  in  the  current  itera- 
tion. Experience  indicates  that  consuming  the  available  space  in  a  side 
constraint  gradually  over  a  moderate  number  of  iterations  is  better  than 
using  it  all  in  one  iteration.  This  gives  the  LP  solver  an  opportunity  to  react 
to  the  rounding  decisions,  and  to  provide  an  improved  basis  for  making 
additional  rounding  decisions. 

When  road  construction  comprises  less  than  10  percent  of  the  side  con- 
straint impact,  "actual  space"  is  adjusted  upward  a  small  amount  to  arrive 
at  "allotted  space."  In  these  cases,  roads  play  a  sufficiently  minor  role  that 
impacts  exceeding  "actual  space"  can  be  accommodated  in  subsequent  LP 
solutions  by  adjustments  in  the  resource  projects  and  traffic  flow  variables 
selected. 

Feasibility  for  Upper  Bound  Constraints — The  "allotted  space"  for  an 
upper  bound  on  constraint  i  iALLOTuB,i)  is  calculated  at  the  beginning  of  an 
iteration  as  follows: 

ALLOTuB.i  =  SPACEi  +  SPACEi  KDELTAiXADJUST)] 

where  SPACEi  is  as  defined  earlier.  DELTAi  and  ADJUST  are  adjustment 
factors  whose  computations  are  described  later  in  this  appendix. 

A  positive  value  for  ALLOTt/g,/  indicates  there  is  space  available  in  side 
constraint  i.  Rounding  options  with  positive  values  for  IMPACT ij  increase 
the  cumulative  total  impact  fi*om  rounding  decisions  (CUMi).  Rounding 
option  j  having  a  positive  IMPACTij  is  within  the  "allotted  space"  if 

IMPACTij  +  CUMi  <  ALLOTuB,i 

where  all  terms  are  as  previously  defined. 

Rounding  options  with  negative  IMPACTij  values  decrease  the  cumulative 
total  impact  fi-om  rounding  decisions  (CUMi).  These  options  are,  therefore, 
always  feasible,  regardless  of  the  value  of  ALLOTubj  or  CUMi. 

Feasibility  for  Lower  Bound  Constraints — ^The  "allotted  space"  for  a 
lower  bound  on  constraint  i  (ALLOTiB.i)  is  calculated  at  the  beginning  of  an 
iteration  as  follows: 

ALLOTLB^i  =  SPACEi  -  SPACEi  [{BELT AiXAD JUST)] 

where  SPACEi  is  as  defined  earlier,  and  DELTAi  and  ADJUST  are  adjust- 
ment factors  whose  computations  are  described  later  in  this  appendix. 

A  negative  value  for  ALLOT i^^.i  indicates  there  is  space  available  in  side 
constraint  i.  Rounding  options  with  negative  values  for  IMPACTij  decrease 
CUMi,  thereby  moving  it  closer  to  the  "actual  space"  for  constraint  i.  Round- 
ing options  with  a  negative  IMPACTij  are  feasible  if 

IMPACTij  +  CUMi  >  ALLOTlb.i 

where  all  terms  are  as  previously  defined. 
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Rounding  options  with  positive  IMPACT ij  values  increase  the  cumulative 
total  impact  CUMi.  They  are,  therefore,  always  feasible,  regardless  of  the 
value  of  ALLOTiB^i  or  CUMi. 

Threshold  Feasibility  for  Equality  Constraints — Equality  con- 
straints have  upper  and  lower  boimds  set  to  the  same  right-hand-side 
value.  Feasibility  with  regard  to  both  boxinds  is  checked.  A  rounding 
decision  is  considered  infe^ible  if  the  calculations  for  either  the  upper  or 
lower  bounds  indicate  an  infeasibiUty. 

Feasibility  for  Constraints  Having  Ranges — Constraints  for  which 
"ranges"  have  been  specified  contain  both  lower  and  upper  bounds.  Feasi- 
bility with  regard  to  both  bounds  is  checked.  A  rounding  decision  is  consid- 
ered infeasible  if  the  calculations  for  either  the  upper  or  lower  bounds 
indicate  an  infeasibility. 

DELTAt — ^This  term  is  an  adjustment  factor  based  on  the  relative  impor- 
tance of  road  construction  in  side  constraint  i.  In  general,  this  factor  makes 
the  threshold  for  constraint  i  more  restrictive  the  greater  the  importance  of 
road  construction  in  constraint  i.  DELTAi  is  calculated  as: 

DELTA.  =  smaller  of  {0.1,  or  -0.5  +  0.05(1/BETA.)} 

BETAi  is  a  measure  of  the  relative  importance  of  road  construction  in  side 
constraint  i,  and  is  calculated  as: 

BETAi  =  Hay  I  /  [Z  \aikXk  \  +1.  laimFm^  +1.  loy  I  ] 
jeS         k  m  jeS 

where 

aij     =    the  coefficient  for  road  project  j  in  constraint  i 

S      =    the  set  of  road  projects  that  have  a  nonzero  activity  level  in 

the  current  LP  solution 
Oik     =    the  coefficient  for  resource  project  k  in  constraint  i 
Xk     =    the  activity  level  for  resource  project  k  in  the  cmrent  LP 

solution 

aim    =    the  coefficient  for  traffic  variable  m  in  constraint  i 
Fm    =    the  activity  level  for  traffic  variable  m  in  the  current  LP 
solution. 

The  range  of  DELTAi  with  respect  to  BETAi  is  shown  below: 


DELTAi 

BETAi 

-0.45 

1.000 

-.43 

.750 

-.40 

.500 

-.30 

.250 

-.17 

.150 

0 

.100 

.10 

<.083 

BETAi,  the  relative  impact  of  road  construction  on  constraint  i,  is  inversely 
related  to  DELTAi.  TMs  means  "allotted  space"  is  a  smaller  percentage  of 
"actual  space"  (SPACE i)  the  larger  the  relative  impact  of  road  construction. 
This  avoids  making  too  many  roimding  decisions  in  one  iteration.  Experi- 
ence indicates  that  better  final  solutions  are  obtained  if  the  number  of 
roimding  decisions  made  in  any  one  iteration  are  somewhat  Hmited.  Limit- 
ing based  on  constraint  i,  however,  becomes  less  important  for  smaller 
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values  for  BETAi .  When  the  importance  of  road  construction  in  constraint  i  is 
less,  the  LP  solution  process  can  more  easily  adjust  to  the  impacts  of  round- 
ing decisions  by  its  selection  of  resource  projects  and  traffic  routing  in  later 
iterations.  In  these  instances,  "allotted  space"  is  less  restrictive  to  avoid 
xmnecessary  iterations  in  the  heuristic  process. 

ADJUST— Tests  indicated  that  the  limits  imposed  by  DELTAi  were  too 
restrictive  in  the  later  iterations.  Too  many  iterations  were  required  to 
achieve  integer  values  for  the  road  projects.  The  factor  ADJUST  was  added 
to  alleviate  this  problem.  This  factor  reduces  the  effect  of  DELTAi  as  the 
number  of  iterations  increase.  ADJUST  is  calculated  as  follows: 

ADJUST  =  larger  of  {0,  or  1.2  -  0.2  iIT-2)} 

where  IT  is  the  number  of  iterations.  IT-2  is  the  number  of  iterations  in 
which  rounding  decisions  have  been  made.  (Iterations  1  and  2  Eire  not 
counted  because  they  do  not  involve  rounding  decisions.)  ADJUST  equals 
zero  beginning  with  iteration  8. 
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